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CHAPTER  I 


STATEMENT  OF  THE  PROBLEM 

No  comprehensive,  test-proven  system  exists  today  for  predicting 
in-soil  tracked  vehicle  performance  quantitatively  and  with  good  accu¬ 
racy.  This  is  not  too  surprising  when  one  considers  the  complexity  of 
the  interaction  of  rotating,  slipping,  geometrically  complex  tracks  of 
variable  tension  and  weight  distribution  with  soils  having  a  wide  va¬ 
riety  of  physical  properties.  The  number  of  terms  required  to  provide 
a  complete  description  of  soil-track  interactions  is  quite  large,  and 
the  task  of  developing  a  theoretically  sound  charaterization  of  all 
major  soil-track  interactions  is  staggering. 

The  need  for  a  quantitative  description  of  the  soil-track  sys¬ 
tem  certainly  exists.  At  present,  the  design  of  tracked  vehicles  is 
based  on  well-defined  knowledge  of  the  operating  characteristics  of  the 
various  components  of  the  vehicle,  but  on  only  vague  knowledge  of  the 
requirements  imposed  on  the  vehicle  by  the  environment  in  which  it  must 
operate,  i.e.  by  off-road  conditions.  Rational  procedure  requires  that 
a  strong  effort  be  made  to  provide  the  tracked-vehicle  designer  with  a 
useful  description  of  the  soil-track  system.  A  first-generation  de¬ 
scription  can  best  be  developed  through  analysis  of  data  obtained  in 
a  laboratory  environment,  where  close  control  and  systematic  program¬ 
ming  of  values  of  the  soil  and  track  variables  are  possible.  A 


1 


2 


photograph  of  the  model  track,  dynamometer  carriage,  and  soil  test  pit 
used  at  the  U.  S.  Ariry  Engineer  Waterways  Experiment  Station  (WES)  in 
the  study  of  soil-track  interactions  is  shown  in  Figure  1.1.  Because 
the  number  of  terms  needed  to  provide  a  full  description  of  all  types 
of  track  behavior  in  all  types  of  soils  is  much  too  large  to  deal  with 
in  a  single  investigation,  an  overall  description  must  be  obtained  in 
parts.  As  a  first  step,  this  thesis  will  develop  a  technique  for  pre¬ 
dicting  the  pull  that  a  track  develops  while  operating  in  a  straight 


FRONT  VIEW  OF  MODEL  TRACK  MOUNTED  IN  LARGE-SCALE 
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FIGURE  1.1 
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RELATION  BETWEEN  PULL  AND  SLIP  FOR  TRACKS 
OPERATING  STRAIGHT-LINE  IN  LEVEL  SECTIONS 
OF  AIR-DRY  DESERT  SAND 


-20  0  20  40  60  60  100 


SLIP,  PERCENT 

FIGURE  1.2 

line  (i.e.  no  maneuvering)  at  near- optimum  slip  (20  percent  slip--see 
Figure  1.2)  in  a  level  test  bed  of  air-dry  desert  sand  (Yuma  sand). 
(Slip  is  defined  on  page  8.)  Hereafter,  pull  at  20  percent  slip  will 
be  designated  P2Q  .  The  study  will  include: 

(i)  The  determination  of  the  soil  and  track  variables  that 
influence  track  pull  most  in  a  desert  sand. 

(ii) -  The  development  of  dimensionless  terms  composed  of  func¬ 
tions  of  the  most  important  independent  soil  and  track  variables,  and 
manipulation  of  these  terms  to  produce  a  single  dimensionless  predic¬ 
tion  term. 

(iii)  Application  of  statistical  techniques  to  determine  the 
functional  relation  between  the  prediction  term  and  dimensionless  pull 
(i.e.  pull/load),  and  to  examine  in  detail  the  adequacy  of  this  relation. 


CHAPTER  II 


DETERMINATION  OF  MOST  IMPORTANT  INDEPENDENT  VARIABLES 

A.  Independent  variables  considered 

Seventeen  independent  (controlled)  variables  were  selected  to 
provide  a  reasonably  comprehensive  description  of  the  soil-track  system 
for  track  pull  performance  at  20  percent  slip  in  air-dry  Yuma  sand. 
These  variables  and  the  ranges  of  values  that  can  be  obtained  for  each 
with  the  WES  test  equipment  are  listed  in  Table  2.1.  Only  16  of 

TABLE  2.1 

INDEPENDENT  SOIL-TRACK  VARIABLES  CONSIDERED  AND 
RANGES  OF  POSSIBLE  VALUES  FOR  EACH 


Variable 


Range  of  Values 


Load,  W 

Soil  strength,  G 

Track-ground  contact  width,  b 

Nominal  track-ground  contact 
length,  i  (on  a  hard  surface) 

Angle  of  approach,  a 

Angle  of  departure,  3 


Near  zero  to  27  kN* 

0.5  to  7»5  MN/m^**  (approximately) 
15.2,  30. 5>  and  6l.0  cm 
6l.O  and  121.9  cm 

5.5  to  18.5  deg  with  forward-end 
bogie  fully  retracted;  21.5  to 
33.0  deg  with  it  fully  extended 

Same  ranges  as  for  q-  (based  on 
positions  of  rear-end  bogie) 
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TABLE  2.1  (Concluded) 


Variable  Range  of  Values 


Road-wheel  diameter 

Road-wheel  spacing  (uniform  with 
at  least  one  inner  road  wheel) 

Pressure  (pneumatic)  in  road 
bogies 

Distribution  of  pressure  in 
bogies 

Location  of  drive  sprocket 

Track-shoe  height 

Track-shoe  thickness 

Track-shoe  spacing 

Index  of  track-belt  tension 

Translational  velocity 

Track-frame  trim  angle,  9 


17.8  cm  only  with  present  system 

For  i  =  6l.O  cm:  20.3  cm  only 
For  i  =  121.9  cm:  20.3,  40.6, 
and  6l.O  cm 

0  to  690  kN/m2  («620  kN/m2  = 
maximum  usable) 

Wide  variety  possible 

Front  or  rear  (can  be  described 
in  terms  of  geometric  location 
relative  to  some  fixed  point  on 
the  track) 

1.3,  2.5,  and  5* **1  cm 

0.32  and  0.64  cm 

3.0  cm  (all  shoes  in);  14.2  cm 
(every  other  shoe  removed) 

p 

Gage  range  =  0  to  20,  700  kN/m  ; 
values  of  1380  to  69OO  kN/m2  have 
been  tested  (see  (v)  on  page  6) 

0  to  0.6  m/sec 

0  to  20  deg  has  been  tested 


*  Ki lone w tons. 

**  Meganewtons  per  cubic  meter. 

these  variables  were  considered  in  this  analysis;  road-wheel  diameter 
was  eliminated  because  it  could  be  tested  at  only  one  value,  17.8  cm. 
Each  of  these  variables  and  the  means  used  to  control  its  value  are 
described  by  Turnape  [  1971]  .  Track  geometry  terms  are  illustrated  in 
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Figure  2.1,  where  schematics  of  the  overall  model  track  and  of  a  sec¬ 
tion  of  WES  band- type  track  are  presented.  Variables  not  shown  in 
Figure  2.1  are  described  below: 

(i)  Index  of  soil  strength,  G.  For  essentially  cohesionless 
soils,  penetration  resistance  gradient  G  characterizes  soil  strength 
as  the  average  slope  of  the  near-linear  curve  of  cone  penetration  re¬ 
sistance  versus  depth  of  penetration.  Cone  penetration  resistance  is 

the  force  per  unit  base  area  required  to  move  a  30-deg  right  circular 

2 

cone  of  3*23-cm  base  area  through  the  soil  normal  to  the  soil  surface 
at  a  rate  of  3-05  cm/sec.  G  is  measured  in  that  depth  range  for  which 
changes  in  soil  strength  noticeably  affect  the  performance  of  a  track 
(taken  as  0  to  15  cm  in  this  study). 

(ii)  Load,  W.  Overall  test  load  W  includes  the  weight  of 
the  model  track  plus  load  applied  at  the  load  axle. 

(iii)  Pressure  (pneumatic)  in  road  bogies.  The  road  bogies, 
taken  together,  support  the  overall  test  load.  Each  road  bogie  in¬ 
cludes  a  pneumatic  cylinder  whose  pressure  is  individually  regulated 
and  monitored. 

(iv)  Distribution  of  pressure  in  road  bogies.  Because  pres¬ 
sure  in  the  road  bogies  is  individually  regulated,  an  unlimited  variety 
of  pressure  distributions  in  the  road  bogies  is  possible. 

vv)  Index  of  track-belt  tension.  This  index  is  obtained  from 
measurement  of  the  pressure  that  a  regulated  hydraulic  jack  exerts  in 
an  outward  direction  on  t::e  track’s  idler  sprocket. 

(vi)  Translational  velocity.  Actual  rate  of  track  horizontal 
advance.  The  travel  ratio  is  the  ratio  of  translational  velocity  to 
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SCHEMATICS  OF  THE  WES  MODEL  TRACK 
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TRACK  SHOE 
HEIGHT 


GROUSER  HEIGHT 
TRACK  TREAD 


TRACK  SHOE 
~  C/i 


SPACING 


BELT 

THICKNESS 


top  view 


S»PE  Y»CW 


TRACK  TREAD  / 

END  VIEW 

b.  SECTION  OF  WES  MODEL  BAND-TYPE  TRACK 


FIGURE  2.1 
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the.  theoretical  rate  of  horizontal  advance  (taken  as  ru)  cos  €  ,  where 
r  is  track  drive  radius,  U)  is  the  angular  velocity  of  the  drive 
sprocket,  and  €  is  the  angle  between  the  bottom  of  the  track  and  a 
horizontal  plane).  Slip  is  defined  as  unity  minus  the  travel  ratio. 

(vii)  Track-frame  trim  angle ,  9.  The  angle  between  the  bottom 

of  the  track  frame  and  the  original  soil  surface. 

B.  Screening  important  variables  by  Placket t-Burman  analysis 

Analysis  of  track  pull  performance  is  greatly  simplified  by  first 
determining  which  of  the  16  independent  soil-track  variables  have  the 
most  influence  on  track  pull,  and  then  restricting  the  description  of 
the  soil-track  system  to  that  provided  by  the  most  important  independent 
variables.  To  identify  the  most  important  variables,  a  screening  de¬ 
sign  of  the  type  developed  by  R.  L.  Plackett  and  J.  P.  Burman  [19^6] 
was  used.  Interest  centered  primarily  in  determining  those  variables 
whose  main  effects  on  track  pull  are  greatest,  and  in  accomplishing 
this  with  a  minimum  of  testing.  These  requirements  suggested  use  of  a 
design  in  which  each  variable  is  tested  at  only  two  values.  The 
Plackett -Burman  designs  drastically  reduce  the  number  of  tests  required 
by  a  2^  factorial  design  at  the  expense  of  a  high  degree  of  fractiona¬ 
tion.  Main  effects  are  not  confounded  with  each  other,  but  are  con¬ 
founded  with  two-factor  and  higher  order  interactions.  Independence  of 
estimated  effects  is  maintained  in  the  design,  and  an  estimate  of  the 
variance  of  the  test  variables  is  provided  by  a  pooled  estimate  of 
variance  of  dunmy  variables. 

Several  Plackett-Burman  designs  are  illustrated  in  Figure  2.2. 

The  size  of  the  matrix  is  indicated  by  the  letter  n  ,  where  n 


SOME  PLACKETT-BURMAN  MATRICES 


TEST* 

no.,  n 

1 

2 

3 

4 

5 

6 

7 

8 


VARIABLE 


4 

4 

4 


B 

4 

4 


C 

4 


4 

4 


4 

4 

4 


(E)** 

4 


4 

4 

4 


<  F) 


(G) 


4 

4 

4 


a.  R  =  4  REAL  VARIABLES;  n  =  R  +  4  =  8  TESTS 


TEST 


VARIABLE 


no.,  n 

A 

B 

C 

D 

E 

F 

G 

H 

m 

U) 

(K) 

i 

4 

4 

- 

4 

4 

4 

- 

- 

4 

- 

2 

4 

- 

4 

4 

4 

- 

- 

- 

4 

- 

4 

3 

- 

4 

4 

4 

- 

- 

- 

4 

- 

4 

4 

4 

4 

4 

4 

- 

- 

4 

- 

4 

4 

- 

5 

4 

4 

aft 

- 

4 

- 

4 

4 

a> 

4 

6 

4 

- 

- 

- 

4 

- 

4 

4 

- 

4 

4 

7 

- 

- 

- 

4 

- 

4 

4 

- 

4 

4 

4 

8 

- 

4 

« 

4 

4 

- 

4 

4 

4 

9 

- 

4 

- 

4 

4 

- 

4 

4 

4 

- 

- 

10 

4 

- 

4 

4 

- 

4 

4 

4 

- 

- 

- 

1  1 

- 

4 

4 

m 

4 

4 

4 

- 

- 

- 

4 

12 

- 

- 

- 

- 

- 

- 

- 

- 

- 

• 

b.  R  =  8  REAL  VARIABLES;  n  =  R +  4  =  12  TESTS 


n  =  12  44-444---4- 

n=  16  4444-4-44--4--- 

n=2C  44--4444-4-4--- 
n=  24  44444-4-44--44- 

n=28  4-4444----4---4 


4  4- 

-444-4-44-4 


c.  FIRST  ROW  OF  PLACKETT-BURMAN  MATRIX  FOR  SEVERAL  VALUES  OF  n 


-  =  LOW  LEVEL 
4  =  HIGH  LEVEL 

*  TESTS  ARE  RUN  IN  RANDOM  ORDER 

**  PARENTHESES  AROUND  A  LETTER  INDICATE  A  DUMMY  VARIABLE.  NO  RE¬ 
STRICTION  IS  PLACED  ON  WHICH  COLUMNS  MAY  BE  CHOSEN  TO  REPRESENT 
DUMMY  VARIABLES. 


FIGURE  2.2 
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identifies  the  number  of  rows,  and  n  -  1  is  the  number  of  columns. 
Plackett-Burman  matrices  have  been  developed  for  variables  tested  at  two 
levels  for  n  =  8,  12,  l6,  ...,  84  ,  88  (92  missing),  96 >  100.  Each  col¬ 
umn  heading  without  parentheses  designates  an  actual  variable  of  the 
test  system.  A  column  heading  with  parentheses  names  a  dummy  (ficti¬ 
tious)  variable;  +  and  -  signs  for  dummy  variables  are  ignored  when 
the  actual  testing  is  done.  Entries  in  each  row  represent  the  values 
of  the  variables  associated  with  one  test.  Each  real  variable  is 
assigned  one  practical  high-level  value  (designated  by  a  +)  and  one 
practical  low-level  value  (indicated  by  a  -)  at  which  it  can  be  main¬ 
tained  nearly  constant  during  testing.  The  sign  is  consistent  along 
any  diagonal  from  right  to  left  beginning  at  any  location  along  either 
the  top  row  or  right-most  column,  except  that  the  last  row  has  all 

minus  signs.  Plackett-Burman  matrices  for  all  values  of  n  have  these 

•  * 

characteristics,  so  that  any  particular  Plackett-Burman  matrix  is  speci¬ 
fied  when  only  its  first  row  is  known.  The  first  rows  of  matrices  for 
several  values  of  n  are  shown  at  the  bottom  of  Figure  2.2. 

Ideally,  a  potentially  complex  system  should  be  examined  at  the 
inception  of  its  study  through  use  of  a  Plackett-Burman  matrix  of  a 
size  just  sufficient  to  treat  all  possibly  important  system  variables. 

If  R  equals  the  number  of  system  (i.e.  not  dummy)  variables,  this  re¬ 
quires  at  least  n  =  R  +  4  tests.  This  results  because  analysis  of 
test  data  is  based  on  a  t-test  of  the  difference  between  mean  values  of 
the  test  response  when  the  variable  of  interest  is  tested  at  its  low- 
level  value  and  when  it  is  tested  at  its  high-level  value.  For  a  given 
level  of  significance,  the  percentage  points  of  t  decrease  very 
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quickly  as  the  number  of  degrees  of  freedom  increase  from  one  to  three, 
and  they  continue  to  decrease  much  more  slowly  thereafter.  Thus,  the 
power  of  the  t-test  (i.e,  the  probability  of  rejecting  the  null  hypoth¬ 
esis)  is  greatly  improved  if  at  least  three  degrees  of  freedom  (i.e. 
three  dummy  variables)  are  present  in  the  test  program,  or  if  no  more 
than  n  -  4  of  the  variables  considered  are  system  parameters. 

Unfortunately,  several  factors  associated  with  the  physical  test¬ 
ing  of  the  model  track  did  not  allow  testing  all  1 6  system  variables 
with  one  Plackett-Burman  design.  During  the  early  stages  of  testing, 

WES  had  available  track  shoes  of  only  one  height,  2.54  cm.  Road -wheel 
spacing  could  be  doubled  for  the  long  track  configuration  (j£  =  121.9  cm), 
since  it  included  seven  road  wheels  and  an  even  number  of  spacings;  but 
spacing  could  not  be  doubled  for  the  short  configuration  (£  =  6l.O  cm) 
because  it  had  four  road  wheels  and  an  odd  number  of  openings.  Also, 
corresponding  patterns  of  road -bogie  pressure  distributions  (other  than 
uniform)  cannot  be  obtained  for  tracks  of  different  lengths.  Finally, 
the  chance  of  obtaining  a  large  t-test  value  for  a  particular  variable 
increases  directly  with  the  size  of  the  difference  between  high-  and  low- 
level  values  at  which  the  variable  is  tested  (if  constant  variance  of 
the  system  variables  is  assumed).  Obvious  interactions  among  several  of 
the  l6  system  variables  (load  and  trim  angle,  for  instance)  restrict  the 
range  of  +  and  -  values  possible  for  the  interacting  variables  if 
they  are  tested  together.  This  restricti/e  influence  is  fairly  great 
for  a  design  that  requires  a  different  combination  of  16  +  and  - 

signs  in  each  of  20  different  tests. 

Ultimately,  the  effects  of  the  1 6  real  soil-track  variables  of 
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Table  2.1  (road-wheel  diameter  was  omitted)  were  studied  in  three  sep¬ 
arate  Plackett-Burman  designs.  Unless  it  was  specified  as  one  of  the 

variables  whose  effect  on  P^q  was  being  studied,  each  variable  had 
its  value  maintained  constant  at  the  level  indicated  in  Table  2.2.  The 


TABLE  2.2 

VALUES  OF  THE  1 6  SAND-TRACK  SYSTEM  VARIABLES  WHEN  THEY  WERE 
HELD  CONSTANT  IN  A  PLACKETT-BURMAN  DESIGN 


Constant  Value 

Variable  Name 

Design 
No.  1 

Design 

No.  2 

Design 
No.  3 

Load,  W 

4450  N*  - 

10,000  N 

Soil  strength,  G 

— 

1.5  MN/nr 

3.3  MN/m3 

Track-ground  contact  width,  b 

-- 

15.2  cm 

15.2  cm 

Track-ground  contact  length,  H 

-- 

121.9  cm 

121.9  cm 

(on  a  hard  surface) 

Angle  of  approach,  a 

22  deg 

-- 

22  deg 

Angle  of  departure,  0 

22  deg 

-- 

22  deg 

Road-wheel  spacing  (uniform  with 

20.3  cm 

20.3  cm 

— 

at  least  one  inner  road  wheel) 

Pressure  (pneumatic)  in  road 

276  kN/m2 

— 

621  kN/m2 

bogies 

Distribution  of  pressure  in  road 

Uniform 

-- 

Uniform 

bogies 

Location  of  drive  sprocket 

Rear 

Rear 

-- 

Track-shoe  height 

2.54  cm 

2.54  cm 

-- 

Track-shoe  thickness 

0.32  cm 

-- 

0.32  cm 

Track-shoe  spacing 

3.0  cm 

■i  •» 

3.0  cm 

Index  of  track-belt  tension 

(all  in) 
4l40  kN/m2 

(all  in) 

6890  kN/m2 

Translational  velocity 

0.6  m/sec 

-- 

0.6  m/sec 

Track-frame  trim  angle,  0 

Unrestrained  Unrestrained 

•"  * 

*  Newtons . 


1 6  system  variables  were  assigned  the  low-  and  high-level  values  listed 
in  Table  2.3  when  they  were  tested  in  three  Plackett-Burman  designs. 


I 
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TABLE  2.3 

LOW-LEVEL  (-)  AND  HIGH-LEVEL  (+)  VALUES  ASSIGNED  TO  THE  16  SAND- 
TRACK  SYSTEM  VARIABLES  IN  THREE  PLACKETT-BURMAN  DESIGNS 


Real 

Variable 

Designation 

System  Variable  Name 

Low-Level 
Value  (-) 

High-Level 
Value  (+) 

Plackett-Burman  Design  No.  1 

A 

Soil  strength,  G 

1.5  MN/m3 

3.0  MN/m3 

B 

Load,  W 

3110  N 

6230  N 

C 

Track  width,  b 

15.2  cm 

6l.O  cm 

D 

Track  length,  JL 

6l.O  cm 

121.9  cm 

Plackett-Burman  Design  No.  2 

A 

Angle  of  approach,  a 

22  deg 

30  deg 

B 

Angle  of  departure,  0 

22  deg 

30  deg 

C 

Translational  velocity 

0.3  m/sec 

0.6  m/sec 

D 

Index  of  track-belt  tension 

1380  kN/m2 

4l40  kN/m2 

E 

Track-shoe  spacing 

3.0  cm 

14.2  cm 

F 

Track-shoe  thickness 

0.32  cm 

0.64  cm 

G 

Minimum  pressure  in  road 

103  kN/m2 

24l  kN/m2 

bogies 

H 

Distribution  of  pressure  in 

Same  for  all 

Uniform  in¬ 

road  bogies 

bogies 

crease,  rear 

to  front 

Plackett-Burman  Design 

No.  3 

A 

Track-shoe  height 

1.3  cm 

5.1  cm 

B 

Road -wheel  spacing 

20.3  cm 

40.6  cm 

C 

Drive -sprocket  location 

Front 

Rear 

D 

Track-frame  trim  angle 

0  deg 

20  deg 

Both  the  first  and  third  Plackett-Burman  designs  were  based  on 
the  matrix  in  part  a  of  Figure  2.2;  the  second  design  was  based  on  the 
matrix  in  part  b  of  Figure  2.2.  Values  obtained  for  track  pull  at 
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20  percent  slip  (P2Q)  in  the  three  Plackett-Burman  designs  are  pre 
sented  in  Table  2.4. 

TABLE  2.4 

VALUES  OF  TRACK  PULL  AT  20  PERCENT  SLIP  OBTAINED 
IN  THREE  PLACKETT-BURMAN  DESIGNS  (NEWTONS) 


Plackett- 

Burman 

Design 

No. 

Test  No. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1 

4777 

2820 

2415 

3083 

1490 

1802 

3430 

1112 

2 

2318 

2349 

2509 

1824 

2064 

1744 

1966 

2251 

1957 

1984 

2028 

2340 

3 

3483 

283 

3796 

3123 

2187 

1693 

826 

3549 

mm  mm 

-- 

— 

— 

The  effect  on  the  test  response  of  each  variable  (both  system  and 
dummy)  is  defined  as 


Effect  on 


A 

y  =  L 

i=l 


[2.B.1] 


where 


Va  = 


the  measured  value  of  the  test  response  (P20)  *n  * 


th 


test 


th 


Si  =  +1  if  the  variable  of  interest  is  +  in  the  i  test,  -1  if 

th 

the  variable  of  interest  is  -  in  the  i  test 
n  =  total  number  of  tests  in  the  design  of  interest 
The  variance  of  the  effects  is  estimated  from  Stowe  and  Mayer  [1966]  by 

m 


1 


(Effect  of  dummy  parameter). 

V(effect  on  y)  =  - - -  [2.B.2] 


m 


1 

! 
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where  m  =  the  number  of  dummy  variables  .  A  value  of  the  t-statistic 

to  test  for  the  significance  of  a  particular  effect  is  given  by 

t  =  where  S.E.  is  the  estimate  of  the  standard  error  of  the 

—  S.E. 

effect  provided  by  the  square  root  of  the  variance  defined  above.  Es¬ 
timates  of  effects,  standard  errors,  and  t- values  (each  based  on  three 
degrees  of  freedom)  are  presented  i  Table  2.5.  Entries  for  the  dummy 
variables  are  included  in  Table  2.5  primarily  because  their  jffects 

TABLE  2.5 

EFFECTS,  STANDARD  ERRORS,  AND  t -VALUES  OBTAINED 
IN  THREE  PLACKETT-BURMAN  DESIGNS 


Plackett-Burman  Design  Number 
12  3 


Var¬ 

ia¬ 

ble 

Effect 

°n  P20 

N  t -Value 

Var¬ 

ia¬ 

ble 

Effect 
°n  P20 

N 

t- Value 

Var¬ 

ia¬ 

ble 

Effect 

0n  P20 

N 

t- Value 

A 

519 

2.334 

A 

-128 

-2.522 

A 

139 

0.384 

B 

1823 

8.197 

B 

11 

0.217 

B 

-877 

-2.421 

C 

980 

4.406 

C 

93 

1.833 

C 

164 

0.453 

D 

461 

2.073 

D 

139 

2.739 

D 

-2241 

-6.185 

(E) 

344 

1.547 

E 

-4 

-0.079 

(E) 

508 

1.402 

(F) 

-23 

-0.103 

F 

-22 

-0.433 

(F) 

231 

O.638 

(0) 

-172 

-0.773 

G 

-385 

-7.586 

(G) 

-287 

-0.792 

H 

-56 

-1.103 

(1) 

-86 

-1.695 

(j) 

-18 

-0.355 

(K) 

-2 

-0.039 

Estimate  of 

standard 

Estimate  of 

standard 

Estimate  of  standard 

error  =  222. 4  N 

error  =  50.75  N 

error  =  362.3  N 

are  useful  in  the  following  analysis. 


Comparison  of  results  within  each  design  can  also  be  made 
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visually  by  using  a  half- normal  plot  of  the  type  proposed  by  Daniel 
[  3 959] •  In  Figure  2.3,  the  absolute  value  of  effect  is  plotted  against 
the  appropriate  percentile  value  of  the  half- normal  distribution  for  all 
system  and  dummy  variables  of  the  three  Plackett-Burman  designs.  Ac¬ 
cording  to  theory,  if  the  effects  are  due  solely  to  a  normally  distrib¬ 
uted  random  error,  they  should  tend  to  fall  along  a  straight  line  when 
the  cumulative  empirical  distribution  is  plotted  in  this  manner.  De¬ 
viation  from  a  straight  line  suggests  significance.  In  design  1,  vari¬ 
ables  A,  B,  and  C  are  singled  out  for  attention;  in  design  2,  variable  G; 
and  in  design  3,  variables  B  and  D.  The  use  of  half-normal  plot  analy¬ 
sis  to  interpret  results  from  designs  that  estimate  variance  of  the 
system  variables  by  means  of  only  three  dummy  variables  is  quite 


HALF-NORMAL  PLOT  OF  EFFECT5  FROM  THREE  PLACKETT-BURMAN  DESIGNS 


FIGURE  2.3 
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perilous;  additional  analysis  of  the  test  results  is  called  for. 

Comparison  of  the  significance  of  effects  obtained  in  different 
Plackett-Burraan  designs  on  the  basis  of  t- values,  is  straightforward 
only  if  the  variance  of  the  system  variables  takes  a  common  value  in 

(S.E.  )2 

'  IT1£IX 9 

each  design.  Comparison  of  the  ratio  - =  51*0  for  designs 

1,  2,  and  3  of  Table  2.5  against  tabled  F  =  27.8  for  a  =  0.05 
and  three  degrees  of  freedom  (from  Pearson  and  Hartley  [1958])  indicates 
at  the  5  percent  significance  level  that  differences  do  exist  among 
variances  in  at  least  two  of  the  designs.  Each  estimate  of  standard 
error  is  valid  only  insofar  as  the  assumption  is  satisfied  that  certain 
higher  order  interactions  are  zero.  Thus,  a  larger  value  of  S.E.  for 
a  given  design  may  indicate  that  interactions  of  system  variables  of 
that  design  exert  more  influences  on  the  test  response  than  do  inter¬ 
actions  of  variables  of  a  design  with  smaller  S.E.  ,  all  other  factors 
being  equal.  It  is  not  necessarily  reasonable,  then,  to  scale  increas¬ 
ing  importance  of  system  variables  from  different  designs  according  to 
increasing  values  of  the  t- statistic. 

If  it  is  hypothesized  that  (a)  variance  of  the  effects,  unaf¬ 
fected  by  interactions  of  the  system  variables,  is  relatively  constant 
for  all  three  designs,  and  (b)  the  smallest  estimate  of  S.E.  from  the 
three  designs  reflects  the  standard  error  not  contaminated  by  interac¬ 
tions,  then  t- values  for  the  system  variables  are  produced  as  shown  in 
Table  2.6: 


TABLE  2.6 


t -VALUES  FOR  THREE  PLACKETT-BURMAN  DESIGNS 
WITH  STANDARD  ERROR  ESTIMATED  AT  50.75  N 


Plackett-Burman  Design  Number 

1  2 

3 

Vari- 

Vari- 

Vari- 

able 

t-Value 

able 

t-Value 

able 

t -Value 

A 

10.227 

A 

-2.522 

A 

2.739 

B 

35-921 

B 

0.217 

B 

-17.281 

C 

19.310 

C 

1.833 

C 

3.232 

D 

9.083 

D 

2.739 

D 

44.158 

E 

-0.079 

F 

-0.433 

G 

-7.586 

H 

-1.103 

At  the  5  percent  significance  level,  the  t-value  must  exceed  3*182  (for 
three  degrees  of  freedom)  to  cause  rejection  of  the  hypothesis  that  a 
given  variable  has  no  effect  on  the  test  responses.  Thus,  in  design  1, 
the  effects  of  A,  B,  C,  and  D  are  indicated  important;  in  design  2, 
variable  G;  and  in  design  3>  variables  B,  C,  and  D. 

Another  means  of  comparing  the  effects  obtained  for  all  system 
and  dummy  variables  of  the  three  designs  involves  the  term 
|  effect  on  P2Q| 

- - - ■?-=; —  X  100  .  For  each  variable  within  a  given  design, 

mean  value  of  P2q 

the  average  of  P^'s  variable  always  at  its  high  level  and 

the  average  of  P2q’s  "^e  variable  always  at  its  low  level  (-) 

are  equidistant  from  the  average  value  of  PgQ  for  all  tests  within 

[effect  on  PgJ 
20 


the  design.  Thus, 


mean  value  of  P, 


X  100  indicates  the  change  in 


percent  from  the  average  value  of  P2Q  caused  by  testing  a  variable 

| effect  on  P^J 

at  its  high  and  low  levels.  Values  of  - = - ■  x  100  are 

mean  value  of  PgQ 

presented  in  Table  2.7: 


TABLE  2.7 
| EFFECT  CN  P  I 

VALUE  0F  A.KKA'V  VALL-V?  Pgp  *  100  F0R  EACH  VARIAF<LE  0F 
THREE  PLACKETT -BURMAN  DESIGNS 


1 

Plackett-Burman  Design  Number 

2 

3 

Vari¬ 

able 

Effect  on  f  n 

AverarG  Valu^?  of 

Vari¬ 

able 

Effect  on  P^ 

20  y  100 

Average  Value  of  P20 

Vari¬ 

able 

Effect  on  P„_ 

- — - - --■-  x  100 

Average  Value  of 

A 

19.8 

A 

6.3 

A 

5.9 

B 

69.7 

B 

0.5 

B 

37.0 

C 

37.5 

C 

4.L 

C 

6.9 

D 

17.6 

D 

6.6 

D 

9'*. 6 

(E) 

13.1 

E 

0.2 

(E) 

21.5 

(F) 

0.9 

F 

1.0 

(E) 

9.8 

(0) 

6.6 

G 

18.2 

(c) 

12.1 

1! 

2.7 

(I) 

'*.1 

(J) 

0.9 

00 

0.1 

Average  value  of  P2Q  =  26l6  N 

Average  value  6f  PgQ  =  2111  N 

Average  value  of  PgQ  =  2368  N 

These  values  indicate  that  the  most  important  system  variables  are  B, 

C,  A,  and  D  from  design  1;  G  from  design  2;  and  D  and  B  from  design  3. 
Taken  together,  analyses  based  on  data  in  Figure  2.3  and  in  Tables  2.6 
and  2.7  indicate  that  the  system  variables  that  influence  P2Q  most  are 
B,  C,  A,  and  D  from  design  1;  G  from  design  2;  and  D  and  B  from 
design  3. 

Findings  in  the  above  paragraph  must  be  interpreted  relative  to 
the  present  state  of  knowledge  of  their  influence  on  track  pull  per¬ 
formance  in  sand.  The  variables  indicated  most  important  are  index  of 
soil  strength  G  ,  load  W  ,  track  width  b  ,  track  length  l  (from 
design  l) ;  minimum  pressure  in  road  bogies  (from  design  2);  road-wheel 
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spacing,  track- frame  trim  angle  0  ,  and  possibly  location  of  the 
drive  sprocket  (from  design  3)«  It  is  common  practice  to  keep  road- 
wheel  spacing  as  small  as  practicable.  Track- frame  trim  angle  0  was 
mechanically  controlled  during  the  tests  of  Plack^tt-Burman  design  3> 
this  can  be  done  in  the  laboratory,  but  not  in  the  field.  Control  over 
0  in  field  operations  is  generally  limited  to  design  of  the  tracked 
vehicle's  weight  distribution,  i.e.  its  at-rest  center  of  gravity  (or 
RCG).  Turnage  [1970]  has  demonstrated  that  locating  the  RCG  from  0  to 
0.3  i  forward  of  the  geometric  center-line  influences  negligi¬ 

bly;  values  of  P2Q  are  reduced  by  moving  the  RCG  toward  the  rear. 

Values  of  P__  were  increased  by  locating  the  drive  sprocket  at  the 
20 

rear  (positive  effect  of  parameter  C,  design  3  in  Table  2.5);  most 
present-day  tracked  vehicles  are  rear- sprocket  driven.  Finally, 
values  of  P20  were  increased  by  decreasing  pressure  in  the  road 
bogies  (negative  effect  of  parameter  G,  design  2).  This  result 
agrees  with  findings  of  practically  all  soil-track  investigators,  in 
that  it  suggests  track  pull  improves  when  load  is  transferred  to  the 
soil  with  a  minimum  of  load  concentration.  Thus,  it  can  be  recommended 
that  (a)  road-wheel  spacing  be  minimized,  (b)  0  be  kept  small  by  lo¬ 

cating  the  track’s  RCG  at  or  forward  of  the  track's  geometric  center, 

(c)  the  track's  drive  sprocket  be  located  at  the  rear,  and  (d)  road 
bogies  be  designed  as  flexible  as  practicable.  Attention  will  be  given 
the  relations  among  P2Q  ,  G  ,  b  ,  ,  and  W  in  the  following 

chapters,  because  these  relations  are  not  well  defined  by  existing  test 
results.  In  all  tests  whose  data  are  examined  in  the  remainder  of  this 
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thesis,  values  of  independent  variables  other  than  G  ,  b  ,  i  ,  and  W 
were  maintained  constant  at  the  following  levels: 


Angle  of  approach,  a 

22  deg 

Angle  of  departure ,  0 

22  deg 

Road-wheel  diameter 

17*8  cm 

Road-wheel  spacing 

20.3  cm 

Pressure  in  road  bogies 

27 6  kN/m2 

Distribution  of  pressure  in  bogies 

Uniform 

Locatic.-  of  drive  sprocket 

Rear 

Track-shoe  height 

2.3b  cm 

Track-shoe  thickness 

0.32  rm 

Track-shoe  spacing 

3.0  cm 

Index  of  track-belt  tension 

6890  kN/m2 

Translational  velocity 

0.6  m/sec 

Track-frame  trim  angle 

Unrestrained 

CHAPTER  III 


DEVELOPMENT  OF  TECHNIQUE  TO  PREDICT  TRACK  PULL 

A.  Determination  of  basic 
dimensionless  terms 

Principles  from  similitude  as  described  by  Murphy  [1950]  are  very 
useful  in  developing  relations  that  permit  reliable  prediction  of  proto¬ 
type  behavior  from  observations  of  model  performance.  The  application 
of  these  principles  to  the  sand-track  system  is  described  briefly  and 
informally  in  the  following  paragraphs;  to  obtain  a  rigorous  descrip¬ 
tion  of  the  principles  of  similitude,  the  reader  is  referred  to 
Murphy  [ 1950] . 

A  basic  part  of  a  similitude  study  is  application  of  dimensional 
analysis  to  the  system  of  concern.  The  first  step  is  to  determine  the 
pertinent  variables.  This  was  done  in  Chapter  II  for  the  system  in¬ 
vestigated  in  this  thesis.  Independent  variables  G  ,  b  ,  l  ,  and 
W  (Table  2.1)  were  indicated  to  be  deserving  of  initial  study,  and  they 
are  referred  to  hereafter  as  the  "basic  variables."  The  next  step  is  to 
determine  the  number  of  dimensionless  and  independent  terms  required  to 
express  a  relation  among  the  variables  of  interest.  These  are  called 
Pi  terms  after  the  Buckingham  Pi  Theorem,  which  states  that  the  number 
of  Pi  terms  (s)  equals  the  total  number  of  variables  involved  (n)  minus 
the  number  of  basic  dimensions  of  the  variables  (b) ,  i.e.  s  =  n  -  b  . 
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For  the  soil- track  system  being  studied,  n  =  5  (G  >  b  ,  X  ,  W,  and 
P^q  ;  the  dependent  variable  P^q  is  included  in  the  total).  Basic 
dimensions  are  units  of  the  variables  in  terms  of  mass,  length,  and 
time  (MLT  units);  these  are  given  in  Table  3*1  for  the  variables  of 
interest.  Only  two  basic  dimensions,  M  and  L  ,  are  included  among 


TABLE  3.1 

BASIC  DIMENSIONS  OF  FIVE  TRACK-SAND  VARIABLES 


Variable 

Basic  Dimensions 
(in  Units  of  Mass-Length-Time) 

0 

OJ 

ft. 

M 

G 

m/l3 

b 

L 

X 

L 

W 

M 

the  five  variables,  so  s  =  n-  b  =  5-  2  =  3»  The  relation  among 
the  variables  of  Table  3*1  can  be  expressed  as  rr^  =  where 

F  indicates  only  that  tt^  (which  includes  the  dependent  variable)  is 
functionally  related  to  tt2  and  rr^  •  Development  of  Pi  terms  from 
the  variables  of  Table  3*1  is  simplified  because  P20  customarily  is 

p 

expressed  in  the  dimensionless  ratio  pull/load,  or  ,  and  purely 

linear  terms  can  be  expressed  as  ratios  of  a  characteristic  track  dimen 
sion.  If  track  length  is  selected  as  the  characteristic  dimension, 

p 

and  rig  =  j  .  The  third  Pi  term  must  include  G  and  W  ,  so 
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that  the  mass  dimension  M  can  be  balanced.  Again  with  track  length  X 

g^3 

as  the  characteristic  track  dimension,  tt,  =  -rr-  .  Thus,  the  functional 

j  w  p 

rPO 

relation  among  Pon  ,  G  ,  b  ,  X  ,  and  W  can  be  expressed  as  -rr- 


B.  Development  of  the  basic - 
variable  prediction  term 


^po  /h  r*  3\ 

The  form  of  the  functional  relation  -jj-  =  Flj  ,  ~y~)  can  be 


determined  only  by  analysis  of  data  from  tests  in  which  values  of  depend - 

P  p 

20  b  G 

ent  Pi  term  -rr-  are  monitored  while  values  of  -r  and  -rr-  are  varied. 


r20 

The  upper  plot  in  Figure  3.1  shows  how  values  of  -rr-  responded  to 

p  W 

no  P  h  l 

changes  in  the  value  of  for  j  =  jj-  .  Data  for  the  two  sizes  of 

tracks  represented  (15.2  cm  X  61.0  cm  and  30.5  cm  X  121.9  cm)  may  be 

no  3 

described  by  one  central  relation,  indicating  that  -rr-  is  related  to 

p  W 

20 

—  in  similar  fashion  for  model  and  prototype.  The  lower  plot  in  Fig- 

^20  b 

ure  3.1  presents  the  versus  relation  for  four  values  of  - 


are 


(l  ’  \  *  £  »  and  g)  *  (Ful1  lines  in  both  plots  of  Figure  3*1 
visual  lines  of  best  fit.) 

To  incorporate  the  influence  of  j  in  the  prediction  term  being 

b  G  i3 

sought,  the  relation  of  j  to  was  needed.  In  the  logarithmic 

b  gX^ 

plot  of  Figure  3*2,  the  inverse  of  j  is  related  to  by  the 


expression 


-i  ■  ■#)' 


2/3 


[3-B.l] 


where  K  is  a  constant  of  proportionality.  Raising  both  sides  of  equa¬ 
tion  [3.B.1]  to  the  3/2  power  yields 
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X3'2  k3/2  y  GX3 

372  *—  ’ 


|  x  (bX)3/2  -  k-2/2  . 


[3.B.2] 


The  value  of  K 


-3/2 


changes  as  a  function  of  -77-  .  Equation  [3.B.2] 


states,  then,  that  for  each  particular  value  of 

^20  G  ( b£ ) 

corresponding  value  of  .  — » — -  is  the 


,  there  is  a 


corresponding  value  of  .  — * — ^ -  is  the  dimensionless  basic- 

variable  prediction  term  indicated  by  relations  among  ,  7  ,  and 

Q  D  W  £r 

Qjlu  r20 

,  to  be  able  to  describe  in  terms  of  the  four  independent 

track-sand  variables  G  ,  b  ,  l  ,  and  W  . 

Dimensional  analysis  has  proved  valuable  in  determining  the  form 

of  a  dimensionless,  basic-variable  prediction  term.  This  technique 

does  not  ascertain  the  functional  relation  between  the  prediction  term 

and  the  dimensionless  performance  term  (or  terms)  it  predicts.  Methods 

of  statistics  will  be  applied  in  the  following  chapter  to  ^define  such  a 

rJh  9)^/2 


relation,  and  to  examine  in  detail  the  adequacy  of  — -  in 

w 


predicting  track  pull  at  20  percent  slip. 


1 

i 


CHAPTER  IV 


CRITICAL  EXAMINATION  OF  BASIC- VARIABLE  PREDICTION  TERM 


A.  Differences  between  prediction  equations  for 
various  levels  of  the  independent  variables 


In  Chapter  III,  it  was  determined  that  there  exists  a  functional 
P20  ,  G(bi)3/2 

T  “d  W  ' 


relation  betv/een 


i.e. 


3a., 

W  1 


A  data  plot  of  the  relation  between 


G(bi) 


3/2 


W 


and  SiMl 


w 


3/2 


w 


[4.A.1] 


is  presented  in 


arithmetic  form  in  Figure  4.1.  (Test  data  used  in  the  analysis  of  Sec¬ 
tion  B  of  Chapter  III  and  in  the  remaining  analyses  of  this  thesis  are 
presented  in  Table  4.1.)  The  relation  of  Figure  4.1  takes  near-linear 
form  when  the  abscissa  scale  is  changed  from  arithmetic  to  logarithmic 
(Figure  4.2).  (In  Figure  4.2  and  hereafter,  log  represents  log1Q.) 

The  equation  of  a  straight  line  in  Figure  4.2  is  of  the  form 

,3/2 


A 

P20  .  ,  G  ( bZ ) ' 

—  =  aQ  4-  &1  log 


[4.A.2] 


where 


a0  ~ 


the  value  of 
intercept) 


A 

P20 

W 


at  log 


G(bl) 


3/2 


W 


=  1.0  (i.e.  the 


an  =  the  slope  of  the  line 
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RELATION  OF  TO 
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KH  DATA  ANALYZED  Hi  THE  DKVKHi  KENT  OK  IHKX»K’TIOW  i  KKK 


TAhLK  U . 1 

r,{hl)3/? 
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AND  IN  THE  LINEAR  REGHKl'niON  OK 


Oil 


LOG 


r,Ol)^ 
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W  ,  N 

1  .  N 
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Values  of  and  a^  were  estimated  from  the  linear  equation  that 

minimized  the  sum  of  squared  deviations  of  the  ordinate  variable  about 
the  regression  line.  The  following  least-squares  equation  resulted: 


A  . 

-|2  =  0.204462  +  0.161470  log-Gtef^ — 


This  equation  is  represented  as  a  continuous  line  in  Figure  4.2. 


[4. A. 3] 


Dashed  lines  in  Figure  4.2  define  95  percent  confidence  limits  for 
*20 

the  value  of  -y-  that  would  be  predicted  by  the  regression  equation  for 
a  particular  observed  value  of  log 


W 


.  From  Davies  [  1961] ,  these 


limits  are  given  by 


Y0  -  t0.02^(n-2)  S 


z(x±  -  x)4 


[4.A.4] 


where 


/\  20 

Yq  =  predicted  value  of  — for  X  =  XQ 

s  =  standard  error  of  estimate  of  Y  on  X 

n  =  number  of  data  points 
3/2 

X.  =  log  -  for  the  i^h  data  point  included  in  the  linear 

1  W 


regression 

Xq  ~  a  future  observation  of  log 


G(b£) 


3/2 


W 


Information  from  the  first  line  of  Table  4.2  was  substituted  in  equation 
[4. A. 4],  and  values  0.068  ,  0.066  ,  0.065,  and  0,067  were  obtained  for 


v\2 


1  (xo  -  x) 
^0.025(11-2)  S  x/1  +  «  +  j(x  .  x) 


for  5^  values  of  0,  1,  2,  and  3, 


respectively.  Thus,  the  separation  between  the  dashed  lines  of  Figure 
4.2  increases  as  the  difference  between  values  of  Xq  and  X  increases, 
but  the  size  of  this  increase  is  practically  imperceptible  for  the  range 


L 


DATA  SEPARATED  BY  LEVELS  OF  b  X  /  ,  Z  ,  AliD 


of  X-  values  considered.  The  magnitude  of  this  separation  indicates 
P20 

that  -rr-  can  be  predicted  with  reasonable  precision  for  the  full  range 

^3/2 


of  values  of  log 


that  might  reasonably  be  encountered. 


The  relation  in  Figure  4.2  is  very  useful.  The  pull  that  a 
track  can  develop  for  a  given  combination  of  values  of  G  ,  W  ,  b  , 
and  H  is  often  important  (consider,  for  example,  towing  another  vehi¬ 
cle  or  pushing  a  soil  mass  with  an  attached  blade  in  front).  Also,  from 
a  large  number  of  field  tests,  the  WES  Soils  Division  [1951]  determined 


that  (O.90)  x 


5 0 

\  w 


attainable  on  level  sand  1  very  closely  matches  the 


maximum  slope  that  a  tracked  vehicle  can  negotiate  in  dry,  loose  sands. 

A 

In  a  linear  regression  of  form  =  a^  +  ,  the  as¬ 

sumption  is  made  that  errors  (e^’s)  are  independent,  have  zero  mean  and 

2 

constant  variance  a  ,  and  follow  a  normal  distribution.  Residuals 

obtained  by  subtracting  predicted  Y^  (i.e.  Y^)  from  measured  Y^  can 

be  considered  to  represent  the  observed  errors  if  the  model  is  correct. 

To  obtain  an  indication  of  whether  the  above  assumptions  are  satisfied, 

the  residuals  can  be  plotted  on  normal  probability  paper.  In  this  plot, 

the  abscissa  represents  the  value  of  the  residual;  the  ordinate  value 

1 

th  m  -  2 

of  the  ni  ordered  residual  is  the  — - —  X  100  percentile  of  the  stand¬ 
ard  normal  distribution,  where  n  is  the  total  number  of  residuals 


and  m  takes  values  1,  2,...,  n  .  Figure  4.3  shows  that  the  data 
closely  approximate  a  straight  line  that  passes  through  (0,  50),  thus 
indicating  a  normal  distribution  centered  on  zero.  (Placement  of  the 
full  line  was  done  by  causing  it  to  pass  through  (0,  50)  and  (s,  84), 
where  s  =  O.O33153  ,  from  Table  4.2.)  The  relation  of  Figure  4.3  does 


NORMAL  PROBABILITY  PLOT  OF  RESIDUALS 
FROM  APPLICATION  OF  EQUATION  [4.A.3] 
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not  contradict  the  assumption  that  the  linear  regression  model  is  correct. 

An  important  concern  of  this  chapter  is  a  validation  of  the  re¬ 
sults  from  the  dimensional  analysis  in  Chapter  III.  One  means  of  evalu- 

P20 

ating  whether  the  influence  on  — —  of  each  variable  included  in 

g(U)3/2 


is  adequately  accounted  for  by  the  form  of  this  term  is  to 


20  GlbAP' 

compare  results  of  linear  regressions  of  — g-  on  log  ^ - 

for  groups  of  test  data  separated  by  levels  of  values  of  b  ,  l  ,  G  , 

and  W  .  Accordingly,  the  test  data  were  segregated  according  to 

(a)  the  six  combinations  of  b  X  i  ,  (b)  six  levels  of  r  ,  and 

(c)  nine  levels  of  W  .  A  summary  of  information  obtained  in  fitting 

P20 

least-squares  regression  equations  of  the  form  =  aQ 
G(bA)^^ 

+  a  log  — — -  to  these  subgroups  of  data  is  given  in  Table  4.2. 

x  W 


The  precision  with  which  a  given  regression  line  could  be  used  to 
P20 

predict  for  each  of  the  data  subgroups  is  indicated  to  be  quite 

good  by  the  small  values  obtained  for  standard  error  of  estimate  s 
(ranging  from  0.0162  to  0.04l8).  Of  more  interest  than  the  variation 
about  ecch  least- squares  line  was  the  determination  of  whether  values  of 
intercept  and  slope  for  each  data  subgroup  differed  significantly  from 
values  for  intercept  and  slope  when  all  data  except  that  of  the  data 
subgroup  of  interest  were  considered.  From  Draper  and  Smith  [  1966] ,  the 

s2Ix^ 

variances  of  intercept  and  slope  can  be  estimated  by  - w  and 

nZ(Xi  -  X)2 

2 

g 

— - ,  respectively,  and  are  designated  here  as  V(a  )  and 

l(\  -  xr 

V(a.)  .  If  small  letters  indicate  values  using  any  one  subgroup  of  data, 


37 


and  capital  letters  indicate  values  using  all  the  remaining  data,  the 
hypotheses  E(aQ)  =  E(A  )  and  E(a1)  =  E(A1)  are  tested  at  the  5  per¬ 
cent  significance  level  by  comparing  tabled  values  of  t^ 

|  ao  -  Aol 

against  values  of  the  computed  t-stabistics  t_  =  — - -  and 

"°  VV(ao}  +  V(V 

I  al  "  ^ll 

tn  =  — ■  - -  ,  respectively.  Denominators  in  these  two  terms 

v<v 

reflect  the  assumption  that  the  variances  of  a^  and  A^  and  of  a^ 
and  A.^  are  unequal;  this  seems  reasonable  from  the  much  smaller  num¬ 
ber  of  data  points  used  to  compute  a^  and  a^  than  to  compute  A^ 
and  A1  ,  respectively.  From  Dixon  and  Massey  [1957]?  an  approxima¬ 
tion  to  the  number  of  degrees  of  freedom  appropriate  for  these  t-tests 
is  defined  as 


df  = 


+  V(AX)  1 

2 

2 

-r  + 

[v(Ax: 

0! 

-  2 


n. 


■«x 

for  t  ,  with  X  =  0  and  1  for  our  considerations. 

Results  of  testing  E(aQ)  =  E(Aq)  and  E(a1)  =  E(A1)  are  pre¬ 
sented  in  Tables  4.3  and  4.4,  respectively.  In  Table  4.3?  examination 
of  the  test  data  caused  the  hypothesis  E(aQ)  =  E(Aq)  to  be  rejected 

in  five  cases:  for  b  X  i  values  15.2  cm  X  6l.O  cm,  15.2  cm  x  121.9 

3 

cm  ,  and  30*5  cm  x  121.9  cm  ,  and  for  G  values  0.90  MN/m  £  G 

<  2  MN/m3  and  2  MN/m3  £  G  £  3  MN/m3  .  Thus,  it  appears  that  for 

G(bj03/? 

log  — v — L - 


W 


20 

W 


=  1  ,  predicted  values  of 
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and  for  the  lowest  strength  sand  test  sections  may  be  different  from 
those  of  the  preponderance  of  the  test  data;  a  similar  statement  might 
be  made  for  the  longest  track,  but  the  test  data  do  not  support  this 
argument  as  well. 

In  Table  4.4  the  hypothesis  E(a^)  =  E(A^)  was  rejected  in  only 

computed  t 

one  case,  for  b  x  i  =  15*2  cm  x  121.9  cm  .  The  ratios  - - — 

i0.975(dft  ) 

computed  tp 

and  t - - - ;  took  values  only  slightly  greater  than  one  for  these 

io.975(dfi;l) 

six  cases,  but  were  only  slightly  less  than  one  for  several  others. 

A  more  detailed  examination  of  values  of  regression  line  results  from 
Table  4.2  is  in  order. 

Values  of  standard  error  of  estimate  s  ,  slope  a«  ,  intercept 
a0 

a^  ,  and  the  ratio  —  from  the  regressions  of  each  of  data  subgroups 

2-7  are  plotted  against  track  width  b  in  Figure  4.4  and  against  track 

length  i  in  Figure  *t.  5-  Corresponding  values  from  data  subgroups  8-13 

and  14-22  are  plotted  against  index  of  soil  strength  G  and  load  W  , 

respectively,  in  Figures  4.6  and  4.7*  Values  of  s  ,  ,  and  a^ 

in  Figures  4.4  and  4.7  appear  to  change  in  random  fashion  as  functions 

of  b  and  W  ,  respectively.  In  Figure  4.6,  the  values  of  a^  ,  a^  , 
a0 

and  —  ,  do  not  show  strong  tendencies  to  depart  from  a  random  pat- 
al 

tern.  Values  of  s  also  appear  not  to  be  related  to  the  value  of  G 
if  results  for  the  regression  of  the  data  subgroup  with  6^G<7  are 
ignored.  The  overall  pattern  of  s  versus  G  in  Figure  4.6  certainly 
does  not  suggest  a  tendency  of  s  to  decrease  as  G  increases.  The 
fact  that  s  for  6^G<7  is  based  on  only  nine  data  points  also 


RELATIONS  OF  S,  aQ,a(,  AND  —  TO  TRACK  LENGTH 

G(bj iy 

FOR  REGRESSION  EQUATION  -“=a0+a,LOG  - — 

FOR  DATA  SUBGROUPS  2-7 


TRACK  WIDTH  b 

O  15.2  CM 
&  30.5  CM 

□  61.0  CM 


NOTE:  S  = 


S  =  STANDARD  ERPOR 
OF  ESTIMATE 

a0  =  VALUE  OF  AT 
LOG  G^b9)3/2=|  o 


LOG  ^ 
a,  =  SLOPE 


FIGURE  4.5 


RELATIONS  OF  S,a0,a,,  AND  ^  TO  LOAD  W 

fCa  G  Cbi 
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w  \ 
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FIGURE  4.7 
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tends  to  support  the  contention  that  the  small  value  o^j  s  for  6  £  G  <  7 


was  the  result  of  chance.  ) 

Only  in  Figure  4.5*  where  the  relations  of  s  ,  a^v.  a^  ,  and 
V.  ~  to  track  length  l  are  presented,  do  there  emerge  distinguishable 

X  -  1  J 

patterns.  (Data  plots  and  least-squares  central  lines  on  which  the  data 

< 

of  Figure  4.5  are  based  are  presented  in  Figure  4.8.)  As  track  length 
increased,  the  values  of  s  generally  increased.  This  tendency  is 
opposite  to  what  might  be  expected  solely  from  a  consideration  of  the 
physical  sizes  of  tracks  represented;  i.e.  as  track  length  increases, 
the  track  might  be  expected  to  become  generally  more  stable  insofar  as 
its  ability  to  maintain  a  consistent  level  of  performance  is  concerned. 
If  the  form  of  ^  truly  reflects  the  influent  of  l  on 

P20 

,  random  patterns  of  s  versus  l  should  emerge.  None  of  the 
values  of  s  are  so  large  as  to  cause  concern  that  the  linear  rela¬ 


tions  between  and  log 

W  W 


for  data  subgroups  2-7  are  not 


reasonably  well  defined;  however,  the  pattern  of  s  versus  l  is 
somewhat  unexpected. 

ao 

In  the  upper  three  plots  of  Figure  4.5,  values  of  an  and  — 

al 

generally  increase  and  those  of  a1  decrease  as  the  value  of  i 

a 

increases.  (The  significance  of  ratio  —  will  be  discussed  in  Sec- 

al 

tion  C  of  Chapter  IV.)  In  all  plots  of  Figure  4.5,  the  relations  of 
the  Y-axis  variable  to  l  are  obscured  by  the  presence  of  three  levels 
of  track  width  b  for  each  value  of  l  .  (A  similar  situation  exists 
in  Figure  4.4.)  To  gain  insight  into  how  the  patterns  of  the  middle 
two  plots  of  Figure  4.5  were  influenced  by  the  b  and  l  main  effects 
and  the  b i  interaction,  values  of  a^  and  a^  were  each  treated 


as  results  from  a  two- factor  factorial  experiment.  Table  4.5  presents 


this  arrangement  for  an  .  In  this  design,  p  =  the  number  of  levels 

U 


TABLE  4.5 

TREATING  a^s  AS  RESPONSES  IN  A  TWO-FACTOR  FACTORIAL  EXPERIMENT 


J =  6l.O  cm 

Xg  =  121.9  cm 

b.^  =  15*2  cm 

0.168611 

0.236609 

B1  =  0.405220 

b2  =  30.5  cm 

0.216304 

0.265057 

B2  =  0.481811 

b^  =  6l.O  cm 

0.225989 

0.214054 

B3  =  0.440043 

■m 

G  =  1.327074 

of  b  =  3  ;  q  =  the  number  of  levels  of  t  =  2  ;  and  n  =  the  num¬ 
ber  of  observations  per  entry  ,  taken  as  one.  For  i  rows  and  j 
columns,  values  involved  in  working  toward  an  analysis  of  variance  are: 


3  B2 


2  2 

Correction  Factor  (CF)  for  the  mean  =  — —  =  =  0.293521 

npq  1x3x2 

2 


SS 


B 


=  £  _  CF  =  (0.405220)  4-  (0^481811)*  +  ( 0.440043)  .  0.293521 

=  0.001471 


i=l 


2  l5 


S3, 


=  S  -  CF  =  1-0^04)  +  (0.716170)  .  0>293521  =  0.00l847 
Lu  np  1x3 


j=l 

3  2 


(BL, ,) 


3S 


BL 


=  ^  S  - n^ - (SSg  +  SSL  +  CF)  =  0.001507 


i=l  j=l 


1 

3 
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The  value  of  >  i»e.  variance  associated  with  the 

values  taken  as  entries  in  Table  4.5,  can  be  estimated  by  pooling  the 

values  of  V(aQ)  for  data  subgroups  2  through  7.  From  Table  4.3,  the 

values  of  V(uq)  for  data  subgroups  2  through  7  range  from  0.0002245 

through  0.0014970,  so  that  an  assumption  of  unequal  variances  of  a^'s 

among  these  data  subgroups  is  indicated.  Since  ANOVA  is  robust  to 

departures  from  constant  variance  (Box  [1954]),  an  assumption  of  equal 

variances  should  not  influence  the  analysis  too  greatly.  A  pooled 

It(n  -  lL  V(a  )  ] 

estimate  for  the  V(aQ)’s  then  is  given  by  — 


*  Dk 

k 

where  k  designates  subgroup  number.  Using  data  for  subgroups  2-7  in 
Table  4.3  yields  a  pooled  estimate  of  0.000621.  The  analysis  of  vari¬ 


ance  is  tabulated  in  Table  4.6.  Values  of 


M3 


B 


MS 


MS 


and 


BL 


res 


MS 


are  1.2 


res 


TABLE  4.6 

ANOVA  FOR  a0's  TREATED  AS  RESPONSES  IN  A 
TWO-FACTOR  FACTORIAL  EXPERIMENT 


Source 

SS 

df 

M3 

Total 

Unknown 

158 

CF 

0.293521 

1 

About  G 

157 

B  main  effect 

0.001471 

2 

0.000735 

L  main  effect 

0.001847 

1 

0.001847 

BL  interaction 

0.001507 

2 

0.000754 

Residuals 

152 

0.000621 

and  1.2,  respectively,  and  FQ  1^2)  =  3.0  .  Thus,  the  hypotheses 

that  b  main  effect  and  b4  interaction  do  not  affect  a^  are  not 


1+9 


rejected  at  the  5  percent  significance  level.  The  calculated  F-value  for 

2  0  001847 

the  hypothesis  HQ  =  a£  -  0  is  .  3.0  .  Because 

=  3.9  >  the  hypothesis  that  l  does  not  affect  aQ  is  not  rejected 
at '''the  5  percent  significance  level.  However,  the  close  proximity  of 

MSj 

the  valucToiV  ^  -  -  to  that  of  FQ  ^5(1  152)  ^oes  cause  some  concern 
regarding  the  influence  of  the  l  main  effect  on  values  of  aQ  . 

The  analysis  of  the  influence  of  b  ,  l  ,  and  bJ£  on  a^ 
follows  a  procedure  parallel  to  that  given  above  for  a^  .  Values 
of  a^  arranged  as  responses  in  a  two-factor  factorial  design  are 
listed  in  Table  4.7  as  follows: 


TABLE  4.7 

TREATING  a/s  AS  RESPONSES  IN  A  TWO-FACTOR 
FACTORIAL  EXPERIMENT 


=  6l.O  cm 

=  121.9  cm 

b^  =  15.2  cm 

0.180383 

0.131734 

B  =  0.312117 

bp  =  30.5  cm 

0.154482 

0.l4l4l6 

Bp  =  0.295898 

b^  =  6l.O  cm 

0.151109 

0.144970 

b3  =  0.296079 

L1  =  0.485974 

Lg  =  0. 418120 

G  =  0.904094 

For  i  rows  and  j  columns : 


CF  =  —  =  0.136231 
npq 


3  B? 


ss 


B 


=  y  -- 

Z,  nq 


CF 


1.000(5^7 


i=l 


2  L2 


SS, 


=  S  -  CF  =  0.000767 

L*  np 

J=1 


-  (SS_  +  SS_  +  CF)  =  0.000520 
D  L 


i=l  j--l 

^residuals  es‘t*mate<*  by  assuming  that  the  V(a^)’s  for  data  sub¬ 
groups  2-7  in  Table  4.4  are  equal,  and  pooling  their  values  to  obtain 
O.OOOI83.  The  results  of  ANOVA,  then,  are  shown  in  Table  4.8.  For  the 


hypothesis  Hq  =  ff,  =  0  ,  Fcalculated  is 


TABLE  4.8 

ANOVA  FOR  a1,s  TREATED  AS  RESPONSES  IN  A  TWO-FACTOR 

FACTORIAL  EXPERIMENT 


Source 

SS 

df 

MS 

Total 

Unknown 

158 

CF 

0.136231 

1 

About  G 

— 

157 

B  main  effect 

O.OOOO87 

2 

0.000044 

L  main  effect 

0.000767 

1 

0.000767 

BL  interaction 

0.000520 

2 

0.000260 

Residuals 

— 

152 

0.000183 

0.000767 

0.000182 


4.2 


and  F0>95(1>152) 


3.9 


Thus,  at  the  5  percent 


J 
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significance  level,  the  test  data  cause  rejection  of  the  hypothesis 

that  i  main  effect  does  not  influence  the  value  of  a.  .  Ratios 
MS  MSrl 

.  —  and  take  values  0.2  and  1.4,  respectively,  each  con- 

KS  Mo 

res  res 

siderably  smaller  than  FQ  ^2  152)  s  3*0  . 

Before  contemplating  an  action  based  on  the  analyses  of  the  ef¬ 
fects  of  b  ,  l9  and  bjfc  on  a^  and  a^  ,  it  is  pxudent  to  examine 

in  more  detail  the  arrangement  of  values  of  Bq  and  a^  in  Tables  4.5 
and  4.7,  respectively.  In  Table  4.5,  by  far  the  smallest  value  of  ag 
was  obtained  for  b  =  15.2  cm,  i  -  6l.O  cm.  The  value  for  b  = 

15.2  cm,  l  =  121.9  cm  was  second  largest  of  those  obtained.  In 
Table  4.7  the  largest  value  of  a^  was  obtained  for  b  =  15.2  cm, 
i  =  61.0  cm,  and  the  smallest  for  b  =  15.2  cm,  l  =  121.9  cm.  Thus, 

the  effects  of  b  ,  i  ,  and  b £  on  values  of  aQ  and  of  a^  would  be 

reduced  considerably  if  data  for  b  =  15.2  cm  were  omitted.  Values  of 
aQ  and  a1  taken  as  responses  in  two-factor  factorial  designs  with 
data  for  b  =  15.2  cm  omitted  are  presented  in  Table  4.9. 

TABLE  4.9 

TREATING  a^s  AND  a^s  AS  RESPONSES  IN  TWO-FACTOR  FACTORIAL 
EXPERIMENTS  WITH  DATA  FOR  b  =  15.2  cm  OMITTED 


For  a. 


b2  =  30.5  cm 
b3  =  6l.O  cm 


■ 

^  =  6l.O  cm 

&2  =  121.9  cm 

0.216304 

0.225989 

0.265057 

0.214054 

L1  =  0.442293 


l2  =  0.479m 


Bi  =  0.481361 
B2  =  0.440043 

G  =  0.921404 


f 


/ 
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TABLE  4.9  (Concluded) 


For  a1 


ll  = 

6l.O  cm 

j&2  =  121.9  cm 

b2  =  30.5  cm 
bg  =  61.0  cm 

0.154482 

0.151109 

0.l4l4l6 

0.144970 

B1  =  O.295898 
B2  =  0.296079 

Li  = 

0.305591 

L2  =  0.286386 

G  =  0.591977 

Pooled  estimates  of  V(a^)  and  V(a^)  were  obtained  by  the 

same  procedure  as  before,  and  the  resulting  analysis  of  variance  is 

M3calc 

shown  in  Table  4.10.  Because  all  ratios  r= -  in  Table  4.10 

Jfcres 


TABLE  4.10 

ANOVA'S  FOR  aQ's  AND  a^s  TREATED  AS  RESPONSES  IN  A 
TWO-FACTOR  FACTORIAL  EXPERIMENT  WITH  DATA  FOR 
b  =  15.2  cm  OMITTED 


Source 

SS 

df 

M3 

MScalc 

MSres 

For 

ao 

Total 

Unknown 

102 

CF 

0.212246 

1 

About  G 

101 

B  main  effect 

0.000427 

1 

0.000427 

0.5 

L  main  effect 

0.000483 

1 

0.000483 

0.6 

BL  interaction 

0.000776 

1 

0.000776 

1.0 

Residuals 


98 


0.000790 


TABLE  4.10  (Concluded) 
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MS  , 
calc 

Source 

SS 

df 

MS 

MS 

res 

For 

ai 

Total 

Unknown 

102 

CF 

0.0876093 

1 

About  G 

101 

B  main  effect 

0.0000003 

1 

0.0000003 

0.0 

L  main  effect 

0.0000921 

1 

0.0000921 

0.5 

BL  interaction 

0.0000117 

1 

0.0000117 

0.1 

Residuals 

-- 

98 

0.0002011 

take  values  much  smaller  than  Fq  =  3*9  >  the  test  data  do  not 

cause  rejection  of  the  hypotheses  that  effects  b  ,  l  ,  and  b l  have  no 
influence  on  a^  and  a^  ,  at  the  5  percent  significance  level,  with  . 
data  for  b  =  15.2  cm  ignored. 

No  definite  answer  is  available  to  the  question  of  why  the 
effects  of  £  on  a^  and  i  on  a^  are  considerably  lessened  when 
data  for  the  15.2-cm-wide  track  are  disregarded.  One  possibility  is 
that  the  narrower  the  track,  the  more  track  pull  is  influenced  by  track- 
sand  interactions  along  the  edge  of  the  track,  as  opposed  to  interactions 
beneath  the  track.  Thus,  15.2  cm  might  be  smaller  than  the  lower  limit 
of  width  that  allows  scaling  of  model  to  prototype.  This  argument  cannot 
be  totally  discounted,  but  it  seems  weal1:  on  the  basis  of  parallel  work 
that  WES  has  done  for  tires  in  sand.  Using  essentially  the  same  approach 
as  followed  in  Chapter  III  of  this  thesis,  Freitag  [1965]  developed  a 
dimensionless  prediction  term  that  can  be  used  to  predict  with  good 

p 

20 

accuracy  the  in- sand  performance  of  tires  whose  hard- surface  con¬ 

tact  widths  range  from  5  to  40  cm  and  more. 


5b 


The  analysis  of  data  in  Tables  4.5  through  4.10  provided  insight 
into  the  causes  for  differences  between  values  of  various  parameters  of 
the  linear  regression  =  aQ  +  &1  log  .  The  overall  implica¬ 
tion  was  that  whatever  changes  in  the  form  of  that  might 

w 

Ppn 

lead  to  improvement  in  the  term's  ability  to  predict  would  involve 

w 

b  and  £  .  Accordingly,  several  dimensionless  forms  of  the  basic- 
variable  prediction  term  were  developed  by  holding  values  of  the  ex¬ 
ponents  of  G  and  W  at  +1  and  -1  ,  respectively  (no  consistent 
pattern  of  effects  of  either  G  or  W  on  either  aQ  or  a^  were 
indicated  in  Figures  4.6  and  4.7)  and  varying  the  values  of  exponents 
of  b  and  £  ,  all  the  while  retaining  the  dimensionless  character  of 


the  overall  prediction  term, 
linear  regressions  of  form 


.Table  4.11  presents  results 


=  aQ  +  ax  log 


Gbxl£x2 


obtained  from 


TABLE  4.11 

PpO  ^ 

C<  MPARISON  OF  VALUES  FR^M  SIMPLE  LINEAR  REGRESSIONS  OF  FORM  «  an  ♦  a.  LOG  °  * 

W  0  1  w 


A 

P20 

W 

X1 

1  ao 

X2 

x  X 

.  tl 

8o  ai 

About  the 
Mean 

KYj  -  Y)2 

Sum  of  Square* 

About 

Regression 

t<Yi  •  ’.«>* 

Due  to 
Regression 

-  Y>2 

Standard  Error 
of  Estimate 

s 

Coefficient 

of 

Correlation 

r 

0 

3.00 

0.208758 

0.107590 

1.087021 

0.593041 

O.49398O 

0.061657 

0.674117 

0.50 

2. SO 

0.164791 

0.139930 

1.087021 

0.394233 

O.692788 

0.050271 

0.798328 

1.00 

2.00 

0.159897 

0.16 1602 

1.087021 

0.228824 

O.858196 

0.038299 

0.888535 

1.20 

1.80 

0.172479 

O.164335 

1.087021 

0.190671 

0.896350 

0.034961 

0.908071 

1.40 

1.60 

0.192374 

0.163294 

1.087021 

0.172862 

0.914159 

0.033288 

0.917408 

1.45 

1.55 

0.198277 

0.162483 

1.087021 

0.171560 

0.915461 

0.033162 

0.917701 

1.50 

1.50 

0.204482 

0.l6l407 

1.087021 

0.171459 

0.915562 

0.033153 

0.917751 

1.55 

1.45 

0.210955 

0.160267 

1.087021 

0.172516 

0.914505 

0.033255 

0.017221 

1.60 

1.40 

0.217658 

0.158887 

1.087021 

0.174680 

0.912341 

0.033463 

0.916135 

l.So 

1.20 

0.246061 

0.151858 

1.087021 

0.193198 

0.893823 

0.035192 

0.906790 

2.00 

1.00 

0.275458 

0.143078 

1.087021 

0.224295 

0.862726 

0.037918 

O.890876 

2.50 

0.50 

O.343668 

0.118444 

1.087021 

0.330237 

0.756784 

0.046010 

O.834386 

3.00 

0 

0.396650 

0.095705 

1.087021 

0.441110 

0.645911 

0.053175 

0.770846 

Values  obtained  for  each  of  r  ,  s  ,  and  SS  about  regression  indicate 

x,  x0 

/yu  ±p  ^ 

the  best  arrangement  of  x.  and  x0  in  - - -  as  values  of  x,  and 

-4  £  W  1 
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x0  both  approach  1.5  •  This  result  is  illustrated  in  Figure  4.9  by  the 


relation  of  s  to  .  Thus,  further  support  is  developed  for  the 

general  conclusion  from  Table  4.10,  i.e.  a  major  part  of  variation  in 

^20  bX  ) 

the  relation  of  to  log  »  * -  is  caused  by  data  for  the  nar¬ 

rowest  track.  From  this  and  the  several  other  analyses  in  Section  A 


of  Chapter  IV,  no  strong  evidence  was  developed  to  indicate  that  sub- 

p 

20 

stantial  improvement  in  ability  to  predict  -rj—  could  be  made  by  chang- 

G(b,)^2 

ing  the  form  of  the  basic-variable  prediction  term  ■  ' -  . 

B .  An  alternate  basic -variable  prediction  term 

— -  *  p 

20 

Another  means  of  diagnosing  whether  the  test  data  suggest 
is  predicted  more  accurately  by  a  term  composed  of  multiplicative 


RELATION  OF  s  TO  x,  FOR  REGRESSION  EQUATIONS 

F?o  GbX|  jLx* 

OF  FORM  =  0.0  +  0.,  LOG  — rrp^- 


:  - 
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FIGURE  4.9 


2.0 

Gb*>Jl** 

W 


5  6 


functions  of  G  ,  b  ,  i  ,  and  W  other  than 


G(bl) 


3/2 


W 


is  to  use  a 


multiple  linear  regression  to  develop  a  prediction  equation  of  a  form 

,3/2 


such  that  the  effectiveness  of 


G(bjfc)' 


W 


and  the  result  from  the  mul¬ 


tiple  linear  regression  can  be  compared.  The  linear  regression  equation 

„/l_  *\  3/2 

-  for  all  158  data  points  is 


that  relates  to  log 


A 

p 

20 

W 


W 


=  0.204482  +  0.161470  log 


G(bA) 


3/2 


W 


[4.B.1] 


This  equation  can  be  expressed  as 
A 

P20 

“p  =  0.204482  +  0.161470  log  G  -  0.161470  log  W 

+  0.242205  log  b  +  0.242205  log  l  .  [4.B.2] 


Equation  [4.B.2]  can  be  transformed  to 

=  0.161470  (log  G  -  0.4271)  -  0.161470  (log  W  -  3.8468) 

+  0.242205  (log  b  -  1.4400)  +  0.242205  (log  j l  -  I.989O)  ,  [4.B.3] 

which  is  of  the  form 

(?  -  y)  =  a^  -  X^)  +  a2(x^  -  %,)  +  a3(X3  -  X^  +  a4(X4  -  X4)  .[4.B.4] 

Equation  [4.B.4]  eliminates  the  intercept  term  aQ  by  transposing  the 
test  data  such  that  the  equation  of  the  line  passes  through  (X^  =  0  , 

X2  =  0)X3  =  0,Xif=0;Y=0). 

P20 

Next,  a  multiple  linear  regression  of  dependent  variable 

w 

and  independent  variables  log  G  ,  log  W  ,  log  b  ,  and  log  l  was 
used  to  develop  an  equation  of  a  form  comparable  to  [4.B.2], 
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The  equation  thus  obtained  is 
A 

p 

20 

-p  =  0.177491  +  0.l4l601  log  G  -  0.169701  log  W 

+  0.248757  log  b  +  0.271217  log  j l  .  [4.B.5] 

Data  pertinent  to  the  multiple  regression  used  to  obtain  equation 
[4.B.5]  are  summarized  in  Table  4.12.  Equation  [4.B.5]  can  be  expressed 
in  the  form  of  [4.B.4]  as 

=  0.l4l601  (log  G  -  0.4271)  -  0.169701  (log  W  -  3.8468) 

+  0.248757  (log  b  -  1.4400)  +  O.271217  (log  l  -  1.9890)  .  [4.B.6] 

The  similarity  of  values  of  '’orresponding  coefficients  in  equations 
[4.B.3]  and  [4.B.6]  is  obvious. 

The  residuals  obtained  from  the  multiple  linear  regression  to 
obtain  equation  [4.B.5]  are  plotted  on  normal  probability  paper  in 
Figure  4.10.  A  corresponding  plot  for  equation  [4.B.1]  appears  in 
Figure  4.3.  The  distribution  of  residuals  in  these  two  figures 
can  be  closely  approximated  by  straight  lines  of  almost  identical 
slope  passing  through  (0,  50),  indicating  very  similar  normal  distribu¬ 
tions  centered  on  zero. 

A  more  definitive  means  of  testing  whether  equation  [4.B.5] 

P 

r20 

describes  the  relation  between  and  log  G  ,  log  W  ,  log  b  ,  and 

log  j l  better  than  does  equation  [4.B.1]  is  to  assess  the  improvement 
in  this  relation  caused  by  estimating  five  coefficients  in  equation 
[4 . B . 5 3  versus  two  in  equation  [4.B.1J.  The  partial  F-test  criterion 
from  Draper  and  Smith  [1966]  can  be  applied  at  the  5  percent 
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MULTIPLE 

LOG  G 

TABLE  14.12 

LINEAR  REGRESSION 

,  LOG  b  ,  LOG  l 

*20 

OF  -fr  ON 
w 

,  LOG  W 

Variable 

Mean 

Standard  Deviation 

Log  b 

1.43997 

0.222934 

Log  G 

0.42708 

0.252037 

Log  W 

3.84680 

0.276024 

Log  j t 

1.98895 

0. 141024 

P2cA 

0.48280 

0.083209 

Correlation  Matrix 


Log  b  Log  G  Log  W  Log  1 


Log  b 

1.000000 

-0.084748 

0.375970 

O.1396H 

0.482651 

Log  G 

-0.084748 

1.000000 

-0.212299 

-0.048155 

0.469801 

Log  W 

0.375970 

-0.212299 

1.000000 

0.283270 

-0.273212 

Log  L 

0.139611 

-0.048155 

0.283270 

1.000000 

0.372595 

V" 

0.482651 

0.469801 

-0.273212 

0.372595 

1.000000 

Variance-Covariance  Matrix 

Log  b 

Log  G 

Log  W 

Log  l 

P20 

W 

Log  b 

0.002263910 

-0.000031590 

-0.000129174 

-0.000285055 

-0.000532886 

Log  G 

-0.000031590 

0.000161152 

0.000000800 

-0.000047450 

-0.000009189 

Log  W 

-0.000129174 

0.000000800 

0.000113233 

0.000022095 

-0.000002682 

Log  /, 

-0.000285055 

-0.000047450 

0.000022095 

0.000116374 

-0.000052148 

P2</' 

-0.000532886 

-0.000009189 

-0.000002682 

-0.000052148 

0.000376011 

/\ 

P20 

W 

=  a0  +  ai  b  +  a2 

Log  G  +  a^  Log  W 

+  Log  l 

Term 

Value 

Standard  Error 

t-Ratio 

ao 

0.177491 

0.047581 

3.73 

"l 

0.248757 

0.012695 

19.60 

°2 

0.l4l601 

0.010641 

13.31 

* 

a3 

-O.1697O] 

0.010788 

-15.73 

1 

"4 

0.271217 

0.019391 

13.99 

Coe*ffici  < 

•nt  of  Multiple 

2 

Determination,  R  = 

0.844279 

Coefficient  of  Multiple  Correlation,  R  0.918847 
Standard  Error  of  Estimate,  r.  -  0,03?83'> 

Humber  of  Data  1’ointr. ,  n  ■  158 
degrees;  of  Freedom,  n  -  p  -  l  153 
Real  dm l  Cum  of  Squares  0,1650 


VALUE 


NORMAL  PROBABILITY  PLOT  OF  RESIDUALS 
FROM  APPLICATION  OF  EQUATION  [4.B.5] 


FIGURE  4.10 


6o 


significance  level  by  comparing  the  quantity 


(SS  from  eq  [4.B.1]  -  SS  from  eq  [4.B.5])  +  3df 
x*6  s _ re  s  _ 

MS  from  eq  [4.B.5] 
res 


against  FQ  ^(3  ^53)  =  2*7  *  ssres  is  est:ijIla‘te<i  as  s  (n  -  P  -  l) 

2 

and  M3  as  s  Based  on  values  from  Tables  4.2  and  4.12,  computed 
F  is  I(Looi922il?6)--^.001°7, 8(1^3)?]  ♦  3  =  2.0  .  ^  the  hypothesis 

that  equations  [4.B.1]  and  [4.B.5]  produce  equally  good  estimates  of 

P20 

the  relation  between  and  log  G  ,  log  W  ,  log  b  ,  and  log  i 

is  not  rejected  at  the  5  percent  significance  level. 

In  summary,  comparisons  of  regression  equations  [4.B.1]  and 

[4.B.5]  and  of  coefficients  that  appear  in  equations  [4.B.3]  and 
[4.B.6]  showed  that:  (a)  values  of  each  of  the  four  coefficients 
of  equation  [4.B.3]  closely  match  those  of  corresponding  coefficients 
in  equation  [4.B.6];  (b)  the  distributions  of  residuals  obtained  from 
application  of  equations  [4.B.1]  and  [4.B.5]  are  very  similar;  and 
(c)  no  significant  difference  was  indicated  in  the  ability  of  equa¬ 


tions  [4.B.1]  and  [4.B.5]  to  predict 


20 

W 


Thus,  the  test  data  do  not 


appear  to  cause  rejection  of  the  hypothesis  that  regression  equation 

P 

[4.B.1]  predicts  w  „  . 

G(b£)^' ^ 

basis  of  this  conclusion,  and  because  -  is  indicated  by  con- 

w  p 

20 

sideration  of  similitude  to  be  able  to  predict  the  —  performance 

of  prototype  tracks  operating  under  conditions  similar  to  those  under 

,3/2 


-n—  as  accurately  as  does  equation  [U.B.5]-  On  the 

w 


which  the  model  tracks  were  tested, 


G(bj e)- 


W 


is  retained  as  the  basic- 


variable  prediction  term  for  tracks  operating  in  sand. 


6i 


C.  Relations  among  P^q  ,  G  ,  b  ,  £  ,  and  W 

for  the  maximum  pull  and  zero  pull  conditions 


F*20  Q(bj£)3/2 

The  dimensionless  form  of  the  terms  — —  and  — * — rf - 

W  W 


has  the 


tendency  to  obscure  the  relations  between  dependent  parameter  P^q  and 
independent  parameters  G  ,  b  ,  £  ,  and  W  .  In  particular,  it  is  of 

interest  to  determine  how  values  of  P2q  change  relative  to  those  of 
load  W  for  a  given  combination  of  track  size  and  index  of  soil 
strength  G  .  The  general  form  of  the  linear  regression  equation  relat- 


to  log 


20  .  . 
— =  a0  +  ai 106 


[4.C.1] 


If  the  description  of  a  particular  track-sand  situation  is  given  by  the 
value  G(bX)^2  ,  denoted  hereafter  by  the  letter  k  ,  then 


ao  +  al  106  W 


P20  =  a0W  +  W  l0g  W  "  ^ai  lo8  W)  W 


[4.C.2] 


To  determine  the  value  of  load  W  that  produces  maximum  P2q  (i.e., 

*  A  * 

P2Q  ) ,  the  derivative  of  P^  with  respect  to  W  is  set  equal  to  0 


yielding 


dP20  + 
dW*  "  a0a 


r  #*/&■,  log  e\"| 

j  log  k  -  log  W  +  W  - jp — jl  =  0  .  [4.C.3] 


1 
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Then, 


*  ao/al 

W  =  10  ^ 


X^  . 

e 


[4.C.4] 


ao/ai  i <■ 

If  10  X  —  is  substituted  for  W  in  equation  [4.C.2], 


/  ao/a  1  k\  r  /  ao/al 

P20  =  (; 10  v  X  f )  X  |a0  +  a^log  k  -  ailog  (10  v  ±  X 


Val 


=  ^10  ^  1  x  ^  x  ^log  . 


[4.C.5] 


Taking  the  ratio  of  values  of  P2^  and  W  yields 


II 

°l 

<  p?l 

(  ar/al  k 

llO  ^  1  X  - 

/\ 

* 

w 

a( 

10 

1  k 
X  e 


=  a-^log  e  =  0.4343a1 


[4.C.6] 


For  W  = 


lo/ai 


4 

dP 

20 

dW* 


2  ^ 
d  P20 

vanishes  and  - 

dW 


a^e  log  e 


thus,  the  second  derivative  test  verifies  a  maximum  value  of  P^_  at 

A  a  /a  A  A  20 

*  O'  lk  *  * 

W  =  10  X  —  .  To  illustrate  the  relation  of  P„~  to  W  ,  con- 

e  20 

sider  an  example  where  b  =  15.2  cm,  t  -  121.9  cm,  G  =  1.0  MN/m  ,  and 

P  3/2 

=  0.204482  to  0.l6l470  log  -  .  Figure  4.11  shows  that  the 

W  A  A  W 

"X*  -X- 

relation  of  P  to  W  is  roughly  parabolic  in  shape,  and  that  for 
do  a 

-X* 

the  conditions  of  this  example,  F^q  is  predicted  to  have  a  maximum 


7^-  =  0.4343a;L  =  0.07013  . 


value  at 


RELATION  OF  TO  W  FOR  A  PARTICULAR  VALUE  OF  G(bfi) 


FIGURE  4.11 
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To  obtain  an  estimate  of  the  precision  of  the  value  of 


A 

* 

P20 


W 


that  produces  maximum  ,  recall  from  equation  [4.C.6]  that  the 


P20 


■yr—  that  corresponds  to  maximum  PQq  is  defined  by  0.4343a^  . 

w* 

Davies  [1961]  defines  100(1  -  2a)  percent  confidence  limits  for  slope 
E(a1)  of  a  simple  linear  regression  analysis  as  a^  +  t 


&  -  x)s 

From  Table  4.2  for  data  subgroup  1  (158  data  points),  a^  =  0.161470  , 

I  £(x ,  7zf 

s  =  0.033153  ,  and  sx  =  J^TTJ -  =  0.472936  .  Then, 


yjL(Xi  -  X) 

QiQ3.3j-23 -  3  0. 0055945  j  and  95  percent  confidence  limits  on 


VI57  X  (0.472936) 
?20 


the  value  of  that  produces  maximum  P  _  are  0.4343(0. 161470 ) 


W 


+  1.96  (0.4343  X  0.0055945),  or  0.07013  +  0.00476  . 

It  is  important,  also,  to  estimate  the  value  of  load  Wj  that 
will  cause  the  track  to  become  immobilized,  i.e.  the  W^.  that  will 
produce  P^q  =  0  .  From  equation  [4 .C .23, 


A 

^0 


=  a0  +  ai  log  k  -  a1  log  Wj-  =  0  , 


or 


A  aQ 

log  WT  =  —  +  log  k  , 
1 
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and 

A  ao/a  1 

Wj.  =  10  x  x  k  .  [4.C.7] 

In  Figure  4.11,  =  io1,26638  x  28.233  kN  =  521  kN  . 

Confidence  limits  on  E(W^)  are  more  readily  obtained  by  dealing 

A  A 

with  the  equation  for  log  WT  than  for  WT  .  [Note  from  equation 

A  1  1  a 

g*q 

[4.C.4]  that  log  (W  for  maximum  P  n)  equals  —  +  log  —  ;  thus, 

av  E-j  0 

.A  .  .A  -L 

log  (W  for  maximum  P2q)  and  log  (W  for  P2q  =  0)  differ  only  in  the  de- 

a0 

nominator  of  the  value  added  to  —  .J  From  Davies  [1961],  Fieller*s 

a. 

^  a0 

Theorem  defines  confidence  limits  for  the  true  value  of  a  ratio  —  . 

al 

where  aQ  and  a^  are  unbiased  estimates,  each  assumed  normally  dis¬ 
tributed,  as  follows: 


a. 


a. 


t2  cov  (a^a^  t 


a?  lal 


2  a, 


t2v(ai) 


[cov^a.^]' 


*VV(aJ  -  -7^  cov  (an,aj  +  ~  V(aJ - —  <V(aJ - v(TT 


O'  a. 


0  1  2  '  1' 
ai 


2 

*1 


t2V(a.) 

*  1 - •  [4-C,8] 

*1 


Values  of  V(sq)  ,  V(a^)  ,  and  cov(aQ,a^)  are  defined  from  Draper  and 

Smith  [1966]  as 


c2£X. 


V(a  )  =  - y  V(aJ  =  - 0  ,  and  cov(an,a, )  = 

O  ..ra/v  vr  \<-  Wv  V 


-o2X 


n£(X.  -  X) 


s(xi  -  X)‘ 


E(Xi  -  X)2 
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The  following  values  are  taken  from  Table  4.2  for  data  subgroup  1  (all 
158  data  points) : 


s2  =  an  estimate  of  o2  =  (0.033153) 2  =  0.00109912 


Z(X.  -  X)2  =  s2  X  (n  -  1)  =  (0.472936)2  X  157  =  35.5869 


SX.  =  nX  =  158(1.72367)  =  272.340 
DC?  =  504.541  ^all  with  X  =  log 


w  J  * 


Substituting  these  values  in  the  expressions  for  V(a q)  ,  V(a^)  , 

and  cov(aQ,a^)  produces  the  following  estimates:  V(a^)  =  O.OOOOQ8627  ; 

V(a^)  =  0.000030 886  ,  and  cov(aQ,a^)  =  -O.OOOO53236  .  With  intercept 

a0  =  0.204482  ,  slope  ^  =  0.161470  ,  t  =  t{a/2>n.2)  -  t(0.025,156) 

=  I.96O  ,  and  the  values  given  above  for  V(aQ)  ,  V(a1)  and  cov(aQ,a1), 

/ao\ 

95  percent  confidence  limits  for  El — J  are  computed  from  equa- 

\al/ 

tion  [4.C.8]  to  be 


=  1.27890  +  0.20483  =  1.07407  to  1.48373  . 

A  aQ 

In  the  equation  log  WT  =  —  +  log  k  ,  the  term  k  equals 

1  a. 

3/2 

G(b/)J/  .  Values  of  track  length  and  width,  b  and  l  ,  respectively, 

are  measured  under  static  conditions  cn  a  flat,  rigid  surface  so  that 
practically  no  variation  is  associated  with  their  measured  (estimated) 
values.  Values  of  index  of  sand  strength  G  do  vary  slightly  over  the 
length  of  a  test  bed.  Standard  deviations  of  log  G  measurements  were 
found  to  be  relatively  constant  for  the  158  test  beds  considered  ( see 

Sn  “  l)s2 

Figure  4.12),  and  a  pooled  estimate  of  form  -  for  q  =  158 

I(n  -  1) 
q 


STANOARO  DEVIATION  S 
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VARIATION  IN  MEASUREMENTS  OF  LOG  G 


0.0646  0.1174 


0.1200 


0.05 


0.04 


0.03 


0.02 


O' - 1 — 

0  6  09  1.0 


3.0  4.0  5.0  6.0  7.0  8.0  9.0  10.0 


INDEX  OF  SAND  STRENGTH  G.  ^ 

m3 


FIGURE  4.12 
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test  beds  produced  the  value  s  =  0.000242236  and  s  =  O.OI56  . 

a0 

The  covariance  of  —  and  G  is  quite  difficult  to  determine. 

al 

For  constant  20  percent  slip  tests  of  the  type  reported  herein,  a 

single  value  of  P2Q  is  developed  within  a  given  test  bed,  so  that 

a0 

there  is  no  possibility  of  a  regression  of  —  versus  log  G  for  a 

a  al 

given  test  bed.  Values  of  —  are  1.4,  1.0,  0.7,  0.9,  1.3,  0*9  for 

al 

G  values  in  subgroups  8-13,  respectively,  of  Table  4.2.  Thus,  no 

a0 

strong  association  is  indicated  between  —  and  G  .  To  be  very  con- 

/_  \  V 


servative 


/ao  \  1  /ao\ 

*0  W  \  1  l-o  \  ™ 

— ,  log  Gl  ,  where  cov  I — ,  log  Gl  is  boun< 

h  /  \ai  / 


+  V(  log  G) 


+  2  cov  { — ,  log  Gl  ,  where  cov  I — ,  log  Gl  is  bounded  above  by 

ru y  v*1 _ \ai  /  nT\~ 

N/V/  X  ^log  =  (°*lo45) (O.OI56)  =  O.OOI63  (with  J V^j 

0.20483  a0  v 

=  •  y  ,  as  indicated  by  the  95  percent  confidence  limits  on  —  ). 

/  a0  \  ^ 
Then,  standard  deviation  of  —  +  log  k  I  is  estimated  as 

_  \  al  / 

^0.01442  =  0.120  ,  and  95  percent  confidence  limits  on  E(  log  Wj) 

are  (~  +  log  k]  +  1.960(0.120).  For  the  example  of  Figure  4.10, 

/a  0  \  ' 

I  ^  +  log  k  1  +  1.960(0.120)  =  (1.26638  +  1.45076)  +  0.23520 
=  2.48194  to  2.95234  .  The  range  covered  by  the  95  percent  confidence 


limits  on  E(log  W^)  is  relatively  small;  however,  when  corresponding 

values  l’or  E(W].)  are  considered  (303*3  kN  and  896. 1  kN,  respectively), 

it  is  obvious  that  only  very  broad  estimates  can  be  made  of  VL.  in 

ac/ai 

terms  of  10  and  k  .  This  should  not  be  too  surprising  when  one 

considers  the  wide  spread  in  X-axis  values  as  opposed  to  Y-axis  values 
in  the  relation  of  Figure  4.2,  on  which  equation  [4.C.7]  is  based.  Load 
W  is  included  in  both  the  x-  and  y-axis  terms  of  Figure  4.2;  it  is  not 


until  W  is  expressed  as  a  function  of  parameters  having  estimated 
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variation  (as  is  done  for  log  in  equation  [4.C.7])>  that  esti¬ 
mates  of  the  value  of  load  that  correspond  to  a  particular  level  of 
P2q  can  be  determined. 

While  determination  of  the  load  that  will  produce  zero  pull  is 
an  important  aspect  in  the  overall  description  of  the  soil-track  system, 
immobilization  due  to  overload  in  a  flat,  sandy  soil  appears  not  to  be 
a  practical  problem  for  a  track  that  has  its  RCG  location  at  its  longi¬ 


tudinal  center  line,  as  was  the  case  for  all  tests  reported  herein. 

Neither  the  immobilization  nor  the  maximum  pull  condition  (developed  at 

P20  \  V 

=  0.0713  +  0.00476)  were  achieved  in  any  test  reported  herein;  this, 


in  spite  of  the  fact  that  the  model  tracks  were  tested  to  values  of 

ground  pressure  up  to  132  kN/m  ,  well  in  excess  of  the  approxi- 

2 

mate  upper  limit  100  kN/m  of  existing  tracked  vehicles. 


The  P2Q  versus  W  relation  illustrated  by  example  in  Fig¬ 
ure  4.11  is  considered  quite  reasonable,  since  a  very  similar  P2Q 
versus  W  relation  substantiated  by  test  data  at  both  the  maximum 


pull  and  zero  pull  conditions  is  reported  by  Turnage  [1971]  for  pneu¬ 


matic  tires  in  sand.  In  a  very  broad  general  sense,  then,  pull  at 
20  percent  slip  for  a  track  with  RCG  at  its  center  line  can  be  con¬ 
sidered  to  behave  in  fashion  similar  to  that  of  a  very  large  pneumatic 
tire,  as  might  be  expected  from  a  consideration  of  the  sand-tire  and 
sand-track  contact  patches. 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 

In  Chapter  I  it  was  pointed  out  that  designers  of  tracked  vehicles 
need  a  quantitative  description  of  the  track-soil  system.  The  complexity 
of  this  system  requires  that  the  initial  stage  of  its  study  be  accom¬ 
plished  in  a  laboratory  environment,  where  close  control  arid  systematic 
programming  of  values  of  the  soil  and  track  variables  are  possible.  It 
was  proposed  that,  in  this  thesis,  data  obtained  from  testing  model 
tracks  in  air-dry  sand  in  the  laboratory  facilities  of  the  U.  S.  Army 
Engineer  Waterways  Experiment  Station  be  analyzed  to  develop  a  means  of 
predicting  the  pull  of  prototype  tracks  operating  at  20  percent  slip 
(i.e.  P^q  or  near-maximum  pull)  on  level,  air-dry  sand. 

In  Chapter  II  the  effects  on  P2Q  of  16  track  and  sand  vari¬ 
ables  selected  to  provide  a  reasonably  comprehensive  description  of  the 
track-sand  system  were  estimated  by  analysis  of  three  Placket t-Burman 
test  designs.  This  type  of  design  requires  a  very  small  number  of  tests 
at  the  expense  of  a  high  degree  of  fractionation,  maintains  independence 
of  the  estimated  effects,  and  estimates  the  error  variance  by  a  pooled 
estimate  of  variance  of  dummy  parameters.  The  influence  of  each  of  the 
16  system  variables  of  the  three  Plackett-Burman  designs  was  judged  on 

the  basis  of  (a)  their  t-statistic  values,  (b)  their  values  of 
[effect  on  P2Q| 

- ■ - = - 7 — •= —  ,  and  (c)  the  deviation  of  their  main  effects  from 

average  value  of  Pon 
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a  straight  line  on  a  half -normal  plot.  Based  on  the  overall  analysis 
from  the  Plackett-Burman  designs,  eight  variables  were  judged  to  affect 
P?Q  significantly:  index  of  sand  strength  G  ,  load  W  ,  track  width 
b  ,  track  length  i  ,  pressure  in  road  bogies,  road-wheel  spacing, 
track-frame  trim  angle  9  ,  and  location  of  drive  sprocket.  Existing 
information  with  respect  to  the  manner  in  which  the  last  four  variables 
influence  track  pull  performance  allows  recommendations  that  (a)  road 
bogies  be  designed  as  flexible  as  practicable,  (b)  road-wheel  spacing 
be  minimized,  (c)  trim  angle  be  kept  small  by  locating  the  at-rest  cen¬ 
ter  of  gravity  (KCG)  at  or  forward  of  the  track's  geometric  center  line, 
and  (d)  the  track  be  rear-sprocket  driven.  Four  variables--G  ,  b  , 
i  ,  and  W--were  selected  as  the  basic  sand-track  variables  to  be  con¬ 
sidered  in  this  study. 


In  Chapter  III  principles  of  similitude  were  used  to  develop 


three  dimensionless,  independent  Pi  terms — ,  j  ,  and  — | - to  ex¬ 

press  the  functional  relation  among  P2q  ,  G  ,  b  ,  X  ,  and  W  . 

Laboratory  test  data  were  analyzed  to  demonstrate  that  the  second  and 

rt  ( h  o.  \ 


third  Pi  terms  can  be  combined  to 


G(bX, 


In  Section  A  of  Chapter  IV,  a  simple  linear  regression  determined 


the  least -squares  fit  of 


20  .  .  G(bx; 

—  to 


as  -T7-  =  0.204482 

w 


+  0.161470  log 


G(bX’ 


five  percent  confidence  limits  for  E 


/P20\ 

y-axis  range  of  v~^~J  t  •  A  plot  of  residuals  on  n< 

'  'predicted 

probability  paper  indicated  a  normal- distribution  centered  on  zero. 


for  158  laboratory  test  data  points.  Ninety- 
imits  for  E^-~p^  covered  a  near-constant 

+  0.065  .  A  plot  of  residuals  on  normal 
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To  test  whether 


G(bX) 


3/2 


W 


'20 


adequately  reflects  the  influence  of 


G  ,  b  ,  X  ,  and  W  on  ,  the  158  data  points  were  separated  into 

subgroups  according  to  (a)  six  combinations  of  b  and  X  ,  (b)  six 

levels  of  G  ,  and  (c)  nine  levels  of  W  .  For  each  data  subgroup,  a 

A  / 

^20  G  ( bX ) 

regression  of  the  form  =  a^  +  log  — ^ ^ - 


was  developed.  The 


coefficients  of  these  equations  were  then  compared  to  and 

A 

p 

from  corresponding  equation  =  An  +  An  log 

W  U  JL 


3/2 


w 


,  based  on  all 


test  data  minus  that  from  the  subgroup  of  interest.  The  hypothesis 
Hq  :  E(aQ)  =  E(Aq)  was  rejected  in  five  cases,  and  HQ  :  E(a1)  =  ECA.^) 
in  one  case.  The  narrowest  track  and  the  lowest  strength  sand  test  sec¬ 
tions  appeared  to  cause  values  of  aQ  different  from  those  obtained  for 
the  preponderance  of  the  test  data. 

A  visual  study  was  made  of  changes  in  values  of  aQ  ,  a.^  ,  and 

s  (standard  error  of  estimate)  for  regressions  based  on  the  21  data  sub¬ 
groups.  For  data  separated  by  levels  of  either  G  or  W  ,  no  system¬ 
atic  patterns  were  noted.  Values  of  aQ  and  a.^  for  data  separated  by 
combinations  of  b  and  X  were  each  analyzed  as  responses  in  a  two- 
factor  factorial  design.  In  these  analyses,  MS  residuals  were  de¬ 
scribed  by  pooled  estimates  of  variances  of  aQ  and  a1  ,  respectively. 
Hypotheses  that  b  ,  0  ,  and  bX  have  no  effect  on  a^  were  not  re¬ 

jected  at  the  5  percent  significance  level.  Corresponding  hypotheses 
for  a1  were  not  rejected  for  b  and  bX  ,  but  was  rejected  for  X  . 
Further  examination  of  values  of  a^  and  a^  in  the  two-factor  fac¬ 
torial  format  showed  that  differences  among  values  of  aQ  and  a^  were 


T 
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considerably  reduced  by  omitting  data  for  the  narrowest  track,  width 
b  =  15.2  cm.  This  raised  some  suspicion  that  b  =  15*2  cm  may  be  near 
the  limit  below  which  relations  for  a  model  track  cannot  be  extrapolated 
to  a  prototype. 

The  above  analyses  developed  no  clear  indication  of  how  exponents 
of  b  and  l  might  be  changed  from  I.50  to  improve  the  ability  of 


*L  x2 
Gb  j R  _ 

W 


to  predict 


20 

W 


prediction  term 


X1  x2 
Gb  _ 

W 


20 


Exponents  x^  and  x^  in  the  dimensionless 


were  varied  from  0  to  3»  and  linear  regres¬ 


sions  of  form  -77-  =  a_  +  a,  log 

W  0  1° 


X1  X2 
Gb  d 

W 


were  obtained.  Values  of 


standard  error  s  indicated  the  optimum  form  of  the  prediction  term 

,3/2 


to  be 


G(bl)- 


W 

In  Section  B  of  Chapter  IV ,  no'  appreciable  difference  was  noted 

P20 

between  two  equations  in  regard  to  their  ability  to  relate  -rr-  to 

w 

A 

P20 

G  ,  W  ,  b  ,  and  l  .  These  equations  were:  =  0.204482 

+  0.161470  log  Gib£^2  /t 

'  A 


based  on  a  simple  linear  regression  of 


20 

W 


on  log 


G(bi) 


3/2 


w 


) 


and 


20 

W 


=  0.177491  +  0.141601  log  G 


-  0.169701  log  W  +  0.248757  log  b  +  0.271217  log  X  (obtained  by  a  mul¬ 


tiple  linear  regression  of 


■20 


w  on  log  G , ,  log  W  ,  log  b  , 
log  X).  This  conclusion  was  based  primarily  on  an  F-test  on  the  re¬ 
duction  in  SS  caused  by  estimating  five  coefficients  in  the  mul- 
tiple  linear  regression  versus  two  in  the  simple  linear  regression. 
These  two  equations  were  also  transformed  to  the  same  general  format, 


' 


1 


L 
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f^20  ^20^  _  —  _  _ 

\"W  W~ /  =  bl^og  G  “  log  G)  +  b2^0g  W  “  ^-og  w)  +  b^(log  b  ”  log  b) 

+  b^Clog  i  -  log  i)  ,  and  values  of  each  of  b^  ,  bg  ,  b^  ,  and  b^ 

from  the  simple  and  multiple  linear  regressions  were  found  to  be  very- 

near  ly  equal.  Thus,  the  comparisons  of  Section  B,  Chapter  IV,  gave  no 

G/b#>3/2 

indication  that  the  exponents  of  G  ,  W  ,  b  ,  and  t  in  — ^ - 

should  be  changed,  or  that  significant  improvements  in  the  ability  to 

p 

predict  20  would  result  from  altering  the  form  of  the  equation 

PA  W  3/2 

«7P  =  0.204482  +  0. 161470  log  . 

w  W 

Of  particular  interest  in  the  description  of  track  pull  in  sand 
are  the  relations  among  P^q  ,  G  ,  b  ,  i  ,  and  W  for  the  maximum 
pull  and  zero  pull  conditions.  In  Section  C  of  Chapter  IV,  an  equation 


of  form  =  aQ  +  log 


was  examined  to  obtain  expressions 


'or  P2q  ,  W  ,  and  for  the  maximum  pull  and  zero  pull  conditions 

in  terms  of  an  ,  a.,  ,  k  ,  and  e  (where  k  =  G(bjfc)^2)  •  Ninety-five 

u  1  p 

20 

percent  confidence  limits  were  obtained  for  (a)  the  value  of  that 

produces  maximum  P^q  >  and  (b)  the  values  of  load  W^.  that  produces 

P20  =  °  * 

Overall,  this  investigation  accomplished  its  primary  goal  of 
determining  a  useful  quantitative  relation  between  track  pull  at 
20  percent  slip  and  four  independent  sand  and  track  variables  that  have 
considerable  influence  on  track  pull.  The  adequacy  of  prediction  equa- 

Ppo  G(bJL)^^ 

tion  -77-  =  0.204482  +  0.l6l470  log  '  -  was  examined  by  several 

w  w 


techniques,  each  of  which  indicated  that  no  significant  improvement  re¬ 
sulted  by  altering  the  form  of  this  equation. 

Future  work  should  include  an  expansion  of  basic-variable  dimin- 

Gfbl)3/2 

sionless  prediction  term  ■  -  to  include  functions  of  other  inde- 

pendent  sand-track  variables  that  influence  track  performance  signifi¬ 
cantly.  Also  the  number  of  dimensionless,  dependent  track  performance 
terms  to  be  predicted  should  be  increased  to  include  measurements  of 
track  sinkage  and  torque  input  required.  Finally,  work  parallel  to  that 
described  in  this  thesis  and  to  that  outlined  in  the  above  sentences  of 
this  paragraph  for  tracks  operating  in  sand  should  be  accomplished  for 
tracks  operating  in  fine-grained  soils  (clays  and  silts). 


■» 
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